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Abstract 

Solitary  waves  in  compressible  media  of  finite  depth  and 
infinite  depth  are  studied.   The  critical  speeds  are  first  ob- 
tained from  the  linearized  equations  and  then  confirmed  by  the 
results  of  the  nonlinear  theory.   Explicit  expressions  for  the 
solitary  waves  are  established  by  a  perturbation  scheme  applied 
to  the  nonlinear  equations. 

The  case  of  a  polytropic  compressible  medium  of  finite 
depth  at  rest  in  the  state  of  equilibrium  is  studied  in  Part  I, 
Solitary  waves  in  compressible  medium  of  infinite  depth  are 
investigated  in  Part  II  and  Part  III.   The  former  concerns  two 
isothermal  layers  at  rest  in  the  state  of  equilibrium  separated 
by  a  contact  surface;  the  latter,  an  isothermal  layer  with  non- 
uniform velocity  distribution  at  equilibrium.   It  is  found  that 
solitary  waves  vanish  at  certain  values  of  characteristic  para- 
meters introduced  in  each  case,  and  especially  no  solitary 
v;ave  solution  exists  for  an  isothermal  layer  of  infinite  depth. 
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Part  I.  Polytropic  Compressible  Media 
1.    Introduction 

The  main  purpose  of  this  v;ork  is  to  extend  Peters  and 
Stoker's  scheme  [1]  to  the  study  of  solitary  waves  in  gravi- 
tating, polytropic  or  isothermal  compressible  media  of  infinite 
or  finite  depth.   These  are  steady  two-dimensional  flows  over 
a  plane  level  bottom  with  a  free  surface  which  may  or  may  not 
be  at  infinity.  The  solitary  waves  are  waves  of  permanent 
type  moving  with  constant  velocity  in  the  horizontal  direction, 
and  the  vertical  displacement  of  the  stream  lines  has  only 
a  single  crest  or  trough  and  tends  to  the  equilibrium  state 
at  infinity.  An  interesting  outline  of  the  history  and  many 
physical  aspects  of  the  solitary  wave  problem  have  been  given 
in  [1].  As  a  supplementary  note  to  the  bibliography  cited 
there  v;e  vjould  like  to  mention  that  the  problem  of  a  solitary 
wave  in  an  incompressible  medium  of  non-uniform  primary  veloc- 
ity distribution  has  been  solved  recently  by  Benjamin  [2]. 
He  used  the  vertical  distance  at  equilibrium  as  one  of  the 
independent  variables  in  place  of  the  stream  function. 
\Jlth   an  approach  different  from  [2]  we  shall  introduce  the 
same  independent  variable  to  modify  our  scheme  in  order  to 
study  solitary  waves  in  compressible  media  with  non-uniform 
velocity  distribution  at  equilibrium. 
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In  Part  I  we  shall  consider  the  simplest  case  of  a  soli- 
tary wave  in  a  one-layer  polytropic  compressible  medium  with 
a  view  to  shov;ing  the  general  approach  to  the  problem  of  this 
kind.  By  polytropic  compressible  medium  we  mean  that  for  the 
compr'essible  medium  there  exists  a  relation  betvjeen  pressure  p 
and  density  p  ,  i.e.  p  =  p^  where  n  >  1  .   For  n  =  1  v/e 
call  the  medium  an  isothermal  compressible  medixim  if  the  equa- 
tion of  state  for  a  perfect  gas  is  used.   The  solitary  v;ave 
problem  is  first  formulated  in  Section  2.   In  Section  3  the 
solution  of  the  linearized  equations  predicts  the  value  of  the 
critical  speed  defined  there,  and  in  Section  4  the  solitary 
wave  near  the  critical  speed  is  investigated  by  the  nonlinear 
theory.   The  results  are  discussed  in  Section  5.   It  is  inter- 
esting to  note  that  the  solitary  viave   solution  does  not  exist 
for  n  =  1  and  n  =  5.15  under  the  present  perturbation 
scheme.   The  former  corresponds  to  an  isothermal  layer  of  infi- 
nite depth  and  the  latter,  a  polytropic  layer  of  finite  depth. 
Furthermore,  for  n  >  5. 15  (1  <  n  <  5«15)   the  speed  of  the 
solitary  wave  is  greater  than  (less  than)  the  speed  /gh 
where  h  is  a  characteristic  length  defined  later,  and  g 
the  gravitational  constant,  and  the  solitary  v;ave  is  one  of 
depression  (of  elevation). 
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In  Part  II  solitary  waves  In  two  isothermal  layers  are 
studied.   These  two  layers  are  separated  by  an  interface, 
the  so-called  contact  surface,  across  which  pressure  and 
velocity  are  continuous  but  density  and  temperature  are  sub- 
ject to  a  discontinuity.   In  principle  it  is  not  difficult 
to  extend  the  method  to  the  cases  of  two  or  more  than  two 
layers  in  each  of  which  n  assumes  different  value;  hov;ever, 
since  more  parameters  must  be  introduced,  the  algebraic  cal- 
culations will  become  prohibitive.   For  the  case  studied  in 
Part  II  the  upper  layer  at  equilibrium  is  extended  to  infin- 
ity and  at  constant  temperature  Tp  ;  and  its  lower  layer 

at  equilibrium  is  of  a  finite  height  h  and  at  temperature 

T 
2 
T,  .  I7e  alivays  assume  that  a  s  7jt=-  >  1  . 

Let  r  =  1  +  [exp  — ^ -  1  ]"  ,   vjhere  Pq  ,  p^  are  the 

Pq/Po 

equilibrium  pressure  and  density  at  the  plane  bottom.   Then 
for  a  given  set  of  a  and  r  there  exist  two  critical  speeds. 
Corresponding  to  each  critical  speed  the  domain  a  >  1  ,  r  >  1 
in  the  a,r-plane  is  divided  into  several  subdomains,  in 
each  of  which  the  solitary  wave  may  be  a  wave  of  elevation 
or  depression,  and  its  speed  may  be  greater  or  less  than  the 
critical  speed.  Along  the  boundary  of  these  subdomains  the 
solitary  v;ave  solution  will  not  exist  under  the  present  scheme. 

In  Part  III  we  shall  investigate  solitary  waves  in  an 
isothermal  compressible  medium  with  non-uniform  velocity 
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distribution  In  the  equilibrium  state.   It  Is  worth-while  to 
mention  that  the  solitary  wave  solution  will  not  exist  as  the 
velocity  becomes  uniform.   This  confirms  what  we  have  shown 
in  Part  I. 

Finally,  vje  would  like  to  remark  that  since  compressibil- 
ity plays  an  Important  role  in  the  study  of  the  atmosphere, 
these  results  obtained,  if  relevant,  may  explain  certain 
geophysical  phenom^ena  related  to  gravitational  waves  in  com- 
pressible media. 


2.  Formulation  of  the  Problem 

VJe  assuine  that  a  body  of  polytropic  or  isothermal  com- 
pressible medium  is  supported  by  a  plane  rigid  bottom  and 
has  a  free  surface  on  vjhich  the  pressure  p  is  zero  and 
there  are  no  geometric  constraints.  A  cross  section  of  the 
medium  in  the  equilibrium  state  is  a  horizontal  strip  of 
finite  or  infinite  depth.   Let  us  assiime  that  a  two-dimensional 
wave  of  permanent  type  v/hich  moves  to  the  left  v;lth  velocity 
c  has  been  created  by  some  disturbance  in  the  medium  initially 
at  rest.  We  choose  a  coordinate  system  moving  with  the  vjave 
such  that  the  x-axls  coincides  with  the  bottom  and  the  y-axls 
passes  through  the  crest  for  a  v;ave  of  elevation  and  the 
trough  for  a  wave  of  depression  (Pig.  1).  As  observed  from 
the  coordinate  system  the  wave  is  stationary  and  the  velocity 
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of  the  medlxim  moving  to  the  right  at  infinity  is  c. 

The  steady  state  equations  governing  the  two-dimensional 
motion  of  a  polytropic  or  isothermal  compressible  medium 
lander  a  gravitational  field  are 
the  equation  of  continuity: 


(2.1)     41^  +--iP^  =  0  , 


oy 


the  equations  of  motion 


du  .  cv  5u     1  5' 


~  du  ,  -v  du     I  dp 

P 


(2.2) 


dv  _   g.  _  1.  ^ 
ax  '  '  "Sy   '    S*  "^  ^ 


'v  Bv   ~  Sv        Id' 
u  -r—  +  V       ~ 


the  specifying  equation: 

(2.3)     P/Po  =  (p/Pq)'"  n  >  1  , 
and  the  equation  of  state: 
f(^,?,T)  =  0  , 

vjhere  u(x,y)  ,   v(x,y)  are  respectively  the  horizontal  and 
the  vertical  velocity  component,   p(x,y)  ,  p(x,y)   and  T(x,y) 
are  respectively  the  pressure,  density,  and  temperature,  and 
g  is  the  gravitional  constant.  In  what  follows  we  shall 
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always  assume  that  the  equation  of  state  takes  the  form 

(2.4)  p  =  R  p  T 

as  the  one  for  a  perfect  gas  where  R  is  the  gas  constant, 
For  n  =  1  the  flov;  is  isothermal,  and  for  an  isentroplc 
flow  of  a  perfect  gas  n  ==  7  =  l.^J-  .  At  the  bottom  and 
on  the  free  Lurface  cwo  conditions  are  imposed: 
At  the  bottom  y  :=  0  : 

(2.5)  V  -  0  ; 

on  the  free  surface  SCxg^yg)  ^   0  , 

(2.6)  u  ^  -1-  V  ^.-  =  0  . 


Let  us  first  consider  the  medium  in  the  equilibrium 
state  moving  v/ith  constant  velocity  u  =  c  .  Hereafter 
v/e  use  the  subscript  <»  to  denote  the  variables  in  the  equi- 


librium  state.   From  (2.1)  -g—  =  0  and  p^  =  Poo^^^  • 


(2.2)  only  gives  p  =  p  (y)  which  is  consistent  with  (2.3) 
Finally  from  (2.2)  to  (2.4)  it  is  obtained  that  for  n  >  1 


1 

~  -   ^^   -     n  h' 
Pn 


n 


(2  7)     2_  =  (1  -  JQzlZ)"-! 


n ;  .  :tocI   9;' 


I  <  a 


:'.■'> 


•x* 


~  n     h^    * 

^0 


and  for     n  =  1    , 


Poo 

~  =  exp(-  g)    , 

Po 


(2.8)  ~  =  exp(-  I)    , 


00  0 


whe 


re  Pq  ,  Pq  and  T„  are  respectively  the  pressure,  density. 


Pq 

and  temperature  at  y  =  0  ,  and  h  =  .   Suppose  that 

SPO 

Pp.  ,  p^  are  always  positive  finite.   Then  it  is  seen  from 

(2.7)  that  for  a  polytropic  medium  of  n  >  1  ,   p   ,  p   and 


T   are  equal  to  zero  at  y  ==  — ^  h  and  the  only  possible 
solution  for  y  >  —V  h  is  p   ,  p   and  T  =  0  .   Prom 

''       n*"i  '  00      ^  00  CO 

(2.8)  we  can  also  observe  that  for  an  isothermal  layer 

P  f   P   "^  0  as  y  •->  CO  .   Therefore,  in  the  equilibrium 


n 


state  a  polytropic  layer  is  of  finite  depth  ■^^  h     while 

an  isothermal  layer  is  of  infinite  depth. 

Now  we  introduce  a  stream  function  ^(x,y)  such  that 


(2.9)     i  =  ?u,  -i—fv. 


0^. 


''T- 
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From  (2.5)  and  (2.6)  It  is  seen  that  both  the  bottom  and  the 
free  surface  are  stream  lines,  i.e.  tp  -   const  along  y  =  0 
and  S(x2,y2)  =  0  .  Me   define  tJj{x,0)   =  0   .   The  constant 
value  of  ^'(Xg^yg)  along  S(x2,y2)  ^  ^    * 

(2.10)    ^(^s'-^S^  ""/         (-pvdx  +  pu  dy)  , 

(X3,0) 

gives  the  mass  flujc  across  any  vertical  plane  from  bottom 
to  the  height  y^  per  unit  breadth.  In  the  equilibrium  state 
u  =  c  ,  for  n  >  1  , 


00 

n 


HTll^ 


^  p  n-±  ^ 


0 

and  for  n  =  1  , 

/oo 
pco^oo^y  =  po^^  • 

0 
Hence  along  S(x2,y2)  =  0  , 

(2.11)  f   =  PQCh  . 

The  totality  of  stream  lines  is  given  implicitly  by 

^(x,y)  =  7  =  const  . 

It  is  assi;imed  that  to  each  value  of  7  such  that  0  <  7  <  p^ch 
there  exists  a  unique  solution  of  y  for  the  above  equation,  i.e. 

(2.12)  y  =  f(x,7)  . 
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The  bar  notation  is  used  hereafter  to  indicate  a  function 
of  x,7  .   From  (2.9)  and  (2.12) 


^  y  i\.        ~ f 

(2.13)         y 


^x-^x^.-=. 


Now  for  any  function  <t'(x,y)  =  0(x,7)  , 


"Sx  ~  "3x   "Sy  "Sx  ~  "5x  ~  P  ^  "sy 
(2.14) 

"SF  =  "53^  "53^  "  P  "^  "57  • 


By  using  f,  u,  p,  and  p  as  dependent  variables,  the  equa- 
tions of  motion  can  be  transformed  to  x,7-plane  with  the 
help  of  (2.13)  to  (2.14).   Prom  (2.2),  (2.3),  (2.13)  and 

(2.14)  we  have,  for  0  <  7  <  pQch,  -  «>  <  x  <  +  eo  , 


u(u^  -pvu^)+vpuu^=-  i(p^  -  p  V  p^)  , 


i.e.      u^  =  -  f^p^  +  f^p^  ; 

(2.15) 

and      u(v^  -"pvu^)  +7*^uv^  =  -g-  -(p  u  p^)  , 


i.e.      u  f^x  ■*-  Vx  =  -  P  S  fy  -  Py  . 
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Together  with  (2.15),  we  have 


P  u  f  ^  =  1 


in  place  of  the  equation  of  continuity,  and  the  specifying 
equation 


£-  =/^Y  n  >  1 


0 


\hi 


to  give  a  relation  between  p  and  p  .  The  boundary  condi- 
tions are : 

f(x,0)  =  0  , 
(2.16) 

p(x,pQch)  =  0  . 

However,  the  above  transformation  is  only  a  fomial  one.   It 
is  seen  from  (2.13)  that  as  p^  — ^  0  ,  f  — >  «>  ,  and  the 
transformation  breaks  down.   In  fact,  the  free  surface  in 
the  x,y-plane  is  the  so-called  vacuum  line  and  the  image 
in  the  x,7-plane  becomes  a  branch  line  [3].  All  the  values 
of  the  variables  at  the  branch  line  must  be  defined  as  the 
limiting  values  of  these  variables  when  y  ]   pQch  . 

In  the  following,  we  introduce  the  dimensionless  vari- 
ables 
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h  '         '^ 2 
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Pr 


'0 
where  h  = .  Then  the  equations  (2.3),  (2.13),  (2.15)> 

PqS 
and  (2.16)  become,  for  0<r)<l,   -eo<^<+oo. 


(2.17) 


puf^  =  1  > 


p  =  Xp   , 


f(l,0)  =  0  ,    p(|,l)  =  0  . 


3.  Linear  Theory.   Critical  Speed 

The  solution  of  (2.17)  for  the  equilibrium  state  u  =  1 
is  found  as  f ollovjs : 

1 
Pq  =  A(l-Ti)  ,  Pq  =  (1-T))^ 
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1  _-  1 
'  ft      ^ 
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(TI.S) 


''ir-n   =  n^ 
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12 
n-1 
^0  =  HTlfl  -  <1-^)  ""  1  '   fo^  n  >  1  , 

=  -  log(l-Ti)  ,   for  n  =  1  , 

V/e  assiAme  that  the  wave  motion  is  a  small  disturbance  super- 
posed on  the  equilibrium  state,  and  let 

1 

p  -  -hii-n)   +  p*  ,     p  =  (i-n)'^  +  p*  , 

n-l 
f  =  ^[1  -  (l-Ti)  "  ]  +  f*  ,   for  n  >  1  , 

=  -  logd-n)  +  f*  ,  for  n  =  1  , 

u  =  1  +  u*  . 

Here  we  may  suppose  that  all  the  starred  quantities  are  uniform- 
ly small  for  0<ti<1,   -oo<^<+co;  however  it  is  un- 
likely that  the  same  would  hold  for  the  derivatives  of  f* 
since  f   has  a  singularity  at  r)  =  1  .  For  the  time  being 
let  us  substitute  these  quantities  in  (2.1?)  and  proceed 
formally  to  neglect  all  the  terms  containing  the  second  order 
products  of  the  starred  quantities  and  their  derivatives.  In 
the  Appendix  v/e  shall  examine  whether  our  linearizing  proce- 
dure is  actually  legitimate.  The  linearized  equations  are 

found  to  be 

1 
ul=   -  (1-n)  ^P*  -  Xf*  , 
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1  1 


f^^  =  -  X(l-Ti)"  f^  -  A(l-ri)   ^  p*  -  p^  , 

(3.1)  _1         1  ^ 

U*  =  -  (1-Tl)    "  p*  -   (1-Tl)^  f*  , 


*=^(l-n)^ 


1  -1 


P*=^U-T1J       p 


* 


> 


subject  to  the  boundary  conditions 

f*(e,o)  =  0  ,         p*(^,l)  =  0  . 
Prom  (3.1)  it  is  found  that 

-  1  i  -1 

i  -2 


(5.2)   +  [(^2  -  i)(i-n)-2]fj  =  0  , 


pj  =  [(1-Ti)'  ^  -  i--(l-T^)-l]-l  [(l-n)'^  f^*  -  AfJ]  . 


Let 


f*  =  F(ti)  Q{'i)    , 


then  from  (3.2) 


Gftfi  +  v^G  =  0  , 


0 


iiiot.j  ■ 


■iu-i)j^ 


Fii 
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i  -1  i  -2 


+    (^2   -  H^(l-n)"^   P  =   v2[(l>Ti)"   ^   .  l^(l.Ti)-l]2   p    , 


(3.3) 

v/here     F     must  satisfy  the   boundary  conditions 
P(0)    =  0    , 


1-  i  i 

liia      (l-Ti)t(l-n)        ""  -  ^]"^[(l-n)"  F„   -  AP]   =  0 


The  solution  for  G(4)  Is  seen  as 

G(U  =  A  cos{v^  +  B)  , 

where  A,B  are  arbitrary  constants.  However,  since  we  are 
only  interested  in  finding  the  value  of  the  critical  speed 
l,     defined  as 

(1) 
&   =  lira  A(v)  , 
V— >-0 

for  simplicity,  the  following  asymptotic  method  is  used  to 
solve  the  equation  for  F  in  (3.3)  for  n  >  1  ,  while  a  gen- 
eral discussion  of  the  solution  will  be  deferred  to  the  Ap- 
pendix.  Let  us  suppose  that,  for  small  values  of  v  , 


(1)   For  a  discussion  of  the  definition  of  critical  speeds, 
c.f.  [1]. 
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?      4 
X   =   £   +   V  A^  +  V  A^  +  • • •  , 


P  =  PQ(r!)  +  v^F^(ti)  ■>,-   v'^F2(ri)  +  ... 


The  equation  governing  ^^i'^  )     is  given  by 

i  -1  i  -2 

[(ni  -  (1-Ti)"   ]Pq^^  +  [-  2i(i-n)'^  +  (1  +^)(i--a)"  ]Fq^ 


+  i(i  -l)(l-ri)-2  p  =  0  , 


subject  to  the  boundary  conditions 

Fq(0)  =  0  , 

1-  i  i 

11m  n^(l-Ti)[ni(l-Tl)   ^  -   1][{1-ti)"  P^  -  £Pq]  =  0  . 

The  solution  of  (3.4)  Is  found  as 

n-1  1 


P^  =  C[l  -  ^(1-Ti)  ^    ]    +  D(l-Tl) 


n-1 


where  C  and  D  are  arbitrary  constants.  By  the  condition 
at  r\  =  1    ,   vie   have  D  =  0  ;  and  if  vie   assume  C  ^   0  ,  i.e. 
the  motion  is  other  than  a  parallel  flow,  then  by  Pq^*^^  ~  ^ 
the  critical  speed  is 

£  =  1    , 

It  is  not  difficult  to  find  the  higher  order  approximations 


i-,-r.M  '  '■■ 
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for  the  solution  of  F  ,  however,  since  we  have  determined 
the  value  of  the  critical  speed  vje  v/111  not  proceed  any  farther. 
For  n  =  1  ,  the  equation  for  F     in  (5.3)  becomes 


(1-X"^)F^^  -  2(1-11)-^  (1-A"^)P^  =  v^d-Tl)"^  d-A"^)^  F  . 


(1)  Suppose  A  ^'  1  .  We  find  that,  by  P(0)  =  0  , 


m^  m. 


'"1        '"o 

p  =  c^[(i-n)  ^  -  (i-Ti)  2]  , 

where  C,   is  an  arbitrary  constant,  and 

m^  =  •![-  1  +  (1  +  ^v2(i.A-^))^/2]  , 


m^  =  |[-  1  -  (1  +  4v2(l-A-l))^/2]  . 


This  solution  is  unbounded  at  t)  =  1  .   Therefore,  in  this 
case  either  linear  theory  fails  or  we  must  set  C,  =  0  and 
P  =  0  . 

(2)  Suppose  A  =  1  ,  V/e  obtain  from  the  equation  for  p^ 
in  (5.2)  that 

(l-Ti)P^  -  F  =  0  , 

and  P  =  0  if  F(0)  =  0  .   Therefore,  we  conclude  that  for 
the  case  of  an  isothermal  layer  of  infinite  depth  either  the 
linear  theory  fails  or  the  solution  must  be  identically  equal 
to  zero. 
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4.  Nonlinear  Theory.   Solitary  Wave  Solution 

Let  us  assume  that  a  solitary  wave  moves  with  a  speed 

such  that  A  =  -^  Is  near  some  positive  value  H    ,   which 

c 

is   to  be  determined  later.      The  equations    (2.17)   can  be  writ- 
ten in  the   form,    for     0<ti<1,      _oo<|<+cxi. 


uf^js    +  u^f^    =    (^-A)pf^   -    ^pf,^   -  p^ 


puf^   =  1 


p  =   -    (i-A)p"  +  2,^^   ,        n  >  I 


Let 


the  above  equations  become 


^a  =  -  \^a   ■*-  ^a^  ' 


^(^^aa  -'-  ^-.'a^   ==  ^^S  "  ^^S  "  ^n 


C+.l) 


puf^  =  1  . 


p  =  -  ep"^  +  ip"  ,    n  >  1 


-.U  [  - 


■^  'CSV'!    C  '    '.'if'     '. 
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together  with  the  boundary  conditions 

f(a,0)  =  0  ,    p(a,l)  =  0  . 

We  assume  that  all  the  dependent  variables  can  be  expanded 
In  integral  powers  of  e  ,  i.e. 


(k.2)  <l>(a,Ti,e)  =  YZ  ^^<^irio,T]) 

k  =  0   ^ 


where  <t)(a,Ti,e)  stands  for  f  ,  u  ,  p  ,  and  p  .   Substitu- 
tion of  (4.2)  in  (4.1)  yields  a  sequence  of  equations  for 
f,  ,  u,  ,  p,  ,  and  p.  .   The  equations  for  the  zero-th 
order  approximation  are,  for  0  <  r\   <  1    ,      -oo<a<+<», 


"Oa   "  "   POa^OT]   '^  ^OaPori    ' 


0  =  -   ^Pnfor,   -  P 


0  Or)        ^On    ' 


(4.3) 


Po^o^on  "  ^  ' 


Pq  =  -«pg    >  n  >  1 


subject  to  the  boundary  conditions 

£^{0,0)    =  0  ,     PQ(a,l)  =  0  . 

We  assume  that  Uq  =  1  and  the  solutions  of  the  above  equa- 
tions are 
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V4A) 


1     ^ 
Pq  =  id-n)    ,  Pq  =    (1-ri)"    ,  for     n  >  1    , 

n-1 
^0  =  H?T^^   -    ^^-""^   ""    ^    '  for     n  >  1    , 


=  -  logd-Ti)    ,  for     n  =  1    , 


Uq    =    1     . 


The  equations   for  the  first  order  approximation  are, 
for     0<Tj<l,      -oo<a<+oo. 


^la  =  -   ^O^Pla  •*•  ^laPOTi 


0  =  P0%   -   ^(Po^lTi   +  Pl%^   -  Pin    ' 

PiVoTi  +  Po^i%  +  PoViTi  =  0  ' 


Pi  =  ^^p^Pq^^"^  -  Pq""  , 


subject  to 

f^(a,0)  =  0  ,     P3^(a,l)  =   0  . 

Elimination  of  u^  ,   f^  ,  and  o^  from  (^.5)  yields 
a  simple  equation 

Plarin  =  ° 
for  p,  .  The  solution  of  p-j^  satisfying  Pj_(a,l)  =  0  is 
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(4.6)     p^^  =  a{(a)(l-n) 

where  a^(a)   is  arbitrary.   By  integration  of  (4.6)  with  re- 
spect to  a     we  have 

^1  ^   a^(a)(l-T])  +  b^(-n)  . 

Since  we  have  assumed  that  the  wave  motion  tends  to  its  equi- 
librium state  at  infinity.  I.e.  a^(a)  -^  0     as  a  — >  -  w  , 
b,(Ti)  =  0  .   Thus 

Pl  "  a^(a)(l-Ti)  , 

where  we  assume  that  a,(a)  ^  0  .   It  is  also  obtained  from 
(4.4)  and  (4.5)  that 

ha   =  (Po-q)"'  (^"'Plna  ^-  ^OaPla^ 

n-1 
=  r^  ai(a)[l  -  i(l-Ti)  ^   ]   ,  for  n  >  1  ; 

and  by  the  equilibrium  condition  at  infinity  we  have 

n-1 


^  -  r^  a^(a)  [1  -  id-T!)  "  ]  ,     for  n  >  1  . 


f 


Since  f,(a,0)  =  0  ,   it  follows  that 

&  =  1    ,  for  n  >  1  , 

which  confirms  the  value  for  the  critical  speed  we  have  ob- 
tained by  the  linearized  equations,  and  also  the  assumption 
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of  positive   I'lniteness  of     £   ,      For     n  =  1    ^   vie  have 

f^(a,i-i)   -  0   . 

This  shovjs  that  for  an  isothermal  layer  of  infinite  depth  the 
medium  alvjays  remains  in  the  state  of  equilibrium,  i.e.  at 
rest,  as  already  discussed  in  the  linear  theory.   The 
solution  for  the  first  order  approximation  is  summarized  as 
follows :  for  n  >  1  , 

1 
p^  =  a^(a)(l-Tl)  ,       p^   =  H^^l^^^  +  l](l-Tl)"  , 

f^  =  a^(a)[l  -  (l-n)  ^   ]    ,  u^  =  -  a3_(a)  -  1  . 

The  equations  for  the  second  order  approximation  are, 

for       0    <    T]    <    1     ,       -co    <    a    <    +eo    , 

"2a  =  -    (P2a%   +  Pla^ln  ^    ''   ^^2a^0n   "^  ^Ic^m^    ' 
(4.8) 

^laa   =    (Po^m   -^  Pl%^    -    (PoS   +  Pl^m   "■  P2%^    "  ^2^ 

Po"2^0Ti  +  PoV2n  +  p2"o^on  "'  Po^i^m  "'-  PiVm  "^  Pi"i%  =  ^  • 

1-n 
1    n  r        ,      n(n-l)  ^2„n-2i 
P2  =  H  PO    fP2  •'  Pi  +  PO  -  —2 PlPo   J 

subject  to  the  boundary  conditions 
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By  elimination  of  \i^   ,   T^   ,     and  p^     from  the  above  equations 
we  obtain 

where 


^-1„     ^     [l  ^l-nr„         n(n~l)     2  n-2i) 
-  PO  Pon^jn  ''O      IPl   -       2 PlPo      ^j, 


gi(cT,n)  =  -  Pi/i^  +  fi^Pi^  , 


(4.10) 


2 


g2(a,n)    =  f^^^   +  u^   +  U3_    , 


g,(a,n)   =  -   Pi^i^   -:-  u^ 


Since  gh(a,Ti)  is  a  known  function  of  a     and  t)  ,  let 


then  the  solution  for  Pp   is,  with  aA(cr)  arbitrary. 


(4.11)    Pgjj  =  a^(a)(l-Ti)  +  j   G2^(a,ri')dTT  . 
It  is  also  obtained  from  (4.8)  and  (4.11)  that 
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^20   =   (Pori^"^^P2aTi  +  ^20   "^  P2a%  "  ^1  ^ 


1  r"^ 


Since  fp  (a,0)  =  0  ,  we  have 

(4.12)    G^(a,0)  +  g^^ia.O)    -  j'      G^(a, n  '  )dTi '  -  g^[a ,0)    =  0 

Making  use  of  {k .K)   and  (4.7),  we  obtain  from  (4.12) 
by  some  involved  but  straight  forward  computation  the  follow- 
ing equation: 

mQa-[''(a)   +  m^aj^(a)a^(a)   +m2aj^(a)   =0   . 


where 


rr,     -  -   2n^  +  13n^  -   15n  -!-  5 
'"O  2(n-l)(2n-l)  ' 


3n  -  1 


ni     _ — 

1  n        ' 


n,     -  3n  -   2 
"^2   ~  2n  -   1    • 


For  reasons  given  in  [1],  we  impose  the  conditions 
aj^(-  <»)  =  aj^'(-  0=)  =  0  ,     aj^(O)  =  0  , 
and  the  solution  for  a-,  (a)  is 

Jmp     p  ,     /~~^p 

a,  (0)  = -~  sech^  ^  a   /-  ~  . 

1         m,       2    y   m^ 
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If  the  results  vje  have  found  can  give  an  accurate  approximation 
to  the  solitary  v;ave  in  a  polytroplc  compressible  medium,  then 
by  returning  to  the  x-7  variables,  \<!e   have 

(4.9) 


p  =-pV2(,.^)^_(,.,)  1^^  .ecY^^h   (A-l)  ]  , 


p  ^  ?,(l-n)^[l  -  (1-A)  ^  (1  -  ^  sech2  _g  ,^  (A-1)  )  ]  , 
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where  n=^,   A=-^,  h=  — -  .      In  terms  of  the  vertical 
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distance,      ^     from  the  bottom  at     x  =  -<»   , 

n-1 
'  '^   o  u  dv  =  o.o.hn   -    fi   -  i2zi  ^' 

Jo 


^  =  f  S.u^dy  =  ^^ch[l   -    (1   -^^)^    ]    . 


3-.  Properties  of  the  Solitary  Wave 


We  have  shown  that 


^i"' = -  ^  =«'=h^  ^v^r^ 


v/here  n  >  1  .   For  n  >  1  ,  both  m,   and  mp  are  positive; 
however,  for  1  <  n  <  nQ  =  5*15  ,     hIq  is  positive;  for  n  =  nQ  , 
iUq     is  equal  to  zero;  and  for  n  >  nQ  ,  m^  is  negative. 
Denote  by  I,  ,  Ip  the  open  intervals   (l,nQ)  and  (nQ,oo) 
respectively.   The  results  are  listed  in  the  following  table: 

A-1     V/ave  type 

+         E 
D 

nip 

The  sign  of  A-1  must  be  such  that  —=■   (A-1)  is  always 

m^ 

positive.   The  last  column  of  the  table  indicates  the  wave 
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type  of  the  solitary  wave  where  we  use  E  for  a  wave  of  ele- 
vation and  D  for  a  wave  of  depression.  As  seen  from  the 
expression  for  f  In  (4.9)  the  wave  type  is  determined  by  the 

mp 
Sign  of  (7\-l)  —  .  For  n  =  1  or  n^  ,  the  solitary  wave 

solution  does  not  exist  since  a, (a)  =  0  for  both  cases. 
It  also  follows  from  the  expression  for  f  that  the  maximum 
deviation  of  the  stream  lines  from  the  horizontal  occurs  at 
the  free  surface. 
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Part  II.  Compressible  Media  of  Infinite  Depth 
with  Two  Isothermal  Layers 
1.  Introduction 

In  the  study  of  the  atmosphere,  if  we  neglect  the  effect 
of  the  earth  rotation  and  curvature,  the  atmosphere  may  be 
regarded  as  a  medium  consisting  of  a  series  of  isothermal 
layei's  extended  to  infinity  over  a  plane  level  surface.   Tliere 
exists  a  linear  relation  between  the  pressure  and  the  density 
In  each  layer  if  the  equation  of  state  for  a  perfect  gas  is 
used.   Each  layer  is  separated  from  the  other  by  a  contact 
surface  across  v/hich  pressure  and  velocity  are  continuous  but 
density  and  temperature  are  subject  to  a  jump.  Let  us  sup- 
pose that  the  medium  is  at  rest  initially  and  a  solitary  vjave 
has  been  created  by  some  disturbance.   We  may  choose  a  moving 
coordinate  system  v;ith  respect  to  v;hich  the  flow  becomes  sta- 
tionary.  Now  v;e  have  a  set  of  equations  as  given  in  Part  I 
govern  the  flovj  in  each  layer  and  tv;o  boundary  conditions  at 
the  bottom  and  the  free  surface.  At  each  Interface  there  are 
tvjo  conditions  of  continuity  to  match  the  solutions  for  each 
layer.  However,  the  difficulty  lies  not  in  finding  the  solu- 
tion but  rather  in  the  interpretation  of  the  results.   This 
will  be  seen  later  even  for  the  two-layer  case  vie   are  going 
to  discuss.   Nevertheless,  the  analysis  presented  here  v;ill 
serve  the  purpose  of  Illustrating  the  method  of  approach. 
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We  shall  consider  the  problem  of  solitary  waves  in  a 
compressible  mediuni  of  two  isothermal  layers  separated  by 
a  contact  surface.   The  temperature  of  the  lower  layer  is 
always  assumed  to  be  less  than  that  of  the  upper  layer.   In 
meteorological  terminology  this  situation  may  correspond  to 
the  so-called  "thermal  inversion."  The  problem  is  formulated 
in  Section  2.   In  Section  3  the  analysis  based  upon  the  line- 
arized equations  is  presented  and  a  solution  in  closed  form 
is  obtained.  The  linear  theory  predicts  two  critical  speeds 
for  a  given  equilibrium  state  . 

In  Section  4  v/e  return  to  the  nonlinear  theory  as  we  did 
in  Part  I .  Tlie  two  critical  speeds  obtained  confirm  those 
due  to  the  linear  theory  and  the  coefficients  of  the  solitary 
wave  solution  indicate  a  complicated  flow  pattern  in  relation 
to  the  parameters  defined.   The  discussion  of  the  results 
will  be  given  in  Section  5'   The  method  employed  in  this  part 
is  quite  similar  to  the  one  of  Part  I;  hoviever,    the  analysis 
given  here  is  self-contained  with  minimum  reference  to  the 
previous  results . 

2.  Formulation  of  the  Problem 

Let  us  assume  that  a  mass  of  compressible  mediiun  consisting 
of  tv;o  layers  fills  up  the  whole  upper  half  space.   The  lower 
layer,  supported  by  a  rigid  plane  bottom,  is  at  temperature  T^ 
and  of  an  equilibrium  height  h  ;  the  upper  layer  at  temperature 
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Tp  separated  from  the  lov/er  layer  by  a  contact  surface  is 
extended  to  Infinity  in  the  state  of  equilibrium.   The  pres- 
sure p  is  assumed  to  be  zero  at  Infinity  and  there  are  no 
geometric  constraints.   A  cross  section  of  the  medium  at  equi- 
libriiim  is  just  the  upper  half  plane  as  shown  in  the  Figure 
2.   We  suppose  that  a  wave  of  permanent  type  moving  to  the  left 
with  constant  velocity  c  has  been  created  in  the  medium 
initially  at  rest.  A  coordinate  system  moving  with  the  wave 
is  chosen  such  that  the  x-axls  coincides  with  the  bottom 
and  the  y-axis  passes  through  the  crest  or  the  trough  of  the 
wave  and  is  positive  upward  (Fig.  2).  As  observed  from  the 
coordinate  system,  the  vjave  is  stationary  and  the  medium  in 
the  state  of  equilibriuim  at  Infinity  moves  to  the  right  with 
constant  velocity  c  . 

The  governing  equations  for  the  lower  layer  are 


o(pu)  .  d(pv)  _  ^ 


(2.1) 


~  Su  ,  ~  c)u     1  ^p 
u-str  +  v-^rr  =  -—  "str 

P 


X     "Sy   "  -v  "Sx  ' 


~  ^v  ,  ~  Sv     _   1  Sp 
"  -53^  •■  ^  -SF  ==  *  ^  -  ^  ^  ^ 


0   '^0 


Pn         P 


Where  p(x,y)   is  the  density,   u(x,y)  ,  v(x,y)  are  the  hori- 
zontal and  vertical  velocity  components,  p(x,y)  is  the  pressure, 
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g  is  the  gravitational  constant,  and  p^  ,   p^  are  the  values 
of  p  and  p  at  y  =  0  in  the  equllibriiwi  state.  The  same 
equations  also  hold  Tor  the  upper  layer. 

V/e  first  consider  the  equilibrium  state  of  the  medium  mov- 
ing ivith  constant  velocity  u,  =  c  ,   where  the  subscript  " 
v;ill  always  denote  the  quantities  in  the  equilibrium  state. 
It  is  seen  from  (2.1)  that  all  the  state  variables  are  functions 
of  y  only.   Suppose  that  Pq  ^  Pn  '  ^^®  values  of  p  and 
p  at  y  =  0  ,  are  given,  and  at  the  interface  y  =  h  p  =  p-,  ^ 


p  =  p,   at  y  =  h~  and  p  =  Po  at  y  =  h   ,  at  infinity 


i\,       i\, 


p  =  p  =  0  .   We  obtain  from  (2.1)  that  for  0  <  y  <  h 


'V. 
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and  for  h  <  y  <  « 


'V 


P,,  -  Po  exp[-  -^  (y-h)]  . 
Pi 

Denote  by  V'(^>y)   the  stream  function  such  that 

pu  -  T//^,.  ,       pv  =  -  ^^  , 

then  the  mass  flux  across  any  vertical  plane  from  y  =  0  to 
y  =  h  per  unit  breadth  is  given  by 


J(x,h)  -  ?^(x,0)  =  J   PcAo^y  =  PjCH^ 
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where     H^  -  M^  [exp(M^  )"-^  -  1  ]h  , 


1    gh 


and  the  mass  fliDc  across  any  vertical  plane  from  y  =  h  to 
y  =  00  per  unit  breadth  is 

p  CO 


^  h 


where     H^  =  M^h 


Pi 


2         gh 

V/e  may  choose 

5^(x,0)   =  0   , 

and  it  follows  that 


^(X,co)  =  PgCH^  . 

Let  the  totality  of  stream  lines  be  given  by 

^(x^y)  =  7 

where   0<y<c.;^  0  <  7  <   PpcH^  •   It  is  assvuned  that  for 
each  7     there  exists  one  and  only  one  stream  line  such  that 

y  =  f(x,7)  . 
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If  X  ,  7  are  chosen  as  Independent  variables,  and  f  ,  p  , 
p  and  u  ,  as  dependent  variables  where  the  bar  notation 
indicates  a  function  of  x  and  7  ,   as  vje  did  in  Part  I^ 
vje  obtain 
for  0  <  7  <  p-,cH,  ,     _  CO  <  X  <  00  ^ 


u  =  -  f  D  -;-  f  p  .  , 
X      7^x    x^7 


"^xx  "■   S'/x  =  -  PSfy  -  P^  , 


P  ufy  =  1  . 


Pi        Pi 

for     PiCHjj^    <  7   <  pgcHg   +  PicHj^    ,  -   co  <  x  <   -a    , 

(2.2) 

U  =  -  F  P   -I-  P  P   , 
X      7  X    X  7  ' 


UP   -:-  U  P  =  -  AgP  -  P   , 
XX    XX      ^7    7  ' 


A  U  P^  =  1  , 


P_  ^  A_  . 

Pi       P 


at  the  bottom,   f(x,0)  =  0  , 
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at  7  =  7-^  =  PgCHg  -:-  p-^cH^  ,    ?{^,7q)    =  0  , 
and  at  the  Interface,  y  =  "^q  "^   P^*^^! 

f(X;7Q)  -  P(>^,7q)  ,  p(-,7o)  =  P(x,7o)  • 

If  we  introduce  the  following  dimensionless  variables 


1  PqCH 


PlC  p^ 


"l 


Pl^  P2  1 


M  =  [exp(M^)"-^  -1]-^,       r  =  l+M 


p^    ^2  ^  0'=^ 


s 


the  equations  (2.2)  become 

for  0<'n<l,    _oo<^<c3, 


U/.  =  -  P/.f^  +  fz-P^  :. 
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(2.3) 


uf^^  •'■  ^e^e  =  -  ^P^  -  Pn  ' 


puf^  =  1  . 


p  =  MAp    ; 


for     l<T]<r,      -co<^<oo. 


U^.    =   -   P,,P^   -I-  P.P„ 


UP..    -:-  U.P.    =   -  PAAF^    -   P, 


^^ 
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AUP„   =  a   , 

n 

P  =  MAA  , 
with  the  boundary  conditions^ 

f(e,0)  =  0  ,  P(4,r)  =  0  , 

f(e,l)  =  P(4.1)  ,  p(e,l)  =  P(4>1)  . 

Here  we  note  that  the  dependent  variables  denoted  by  lower 
case  letters  alv/ays  correspond  to  the  lower  layer  and  those 
denoted  by  capital  letters  always  correspond  to  the  upper  layer. 
It  is  also  understood  that  the  values  of  the  variables  at  the 
branch  line  ri  =  r  are  defined  as  the  limits  of  these  variables 
as  T]  I  r  . 
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3.  Linear  Theory.   Critical  Speeds 

For  the  equilibrium  state  of  a  steady  flow  with  u  =  U  ^  1 
we  find  from  (2-5)  tliat 

Pq  =  A{r-Ti)  ,  Pq  =  A(r-Ti)  , 

Pq  =  (M)"^  (r-Ti)  ,  Aq  =  (M)'^  (r-Ti)  , 

fQ  =  -  M  log  ^  ,        Fq  =  -  Ma  log  ^  -  M  log 


r-1 


Let  us  consider  a  small  disturbance  superposed  on  the  equilib- 
rium state  and  we  may  v/rite 

p  =  ?v(r-Ti)  -t-  p*  ,         P  =  A(r-r|)  -:-  P*  , 


p  =  (M)"-'-  (r-r|)  +  p*  ,    A  =  (M)""'-  (r-r))  -i-  A*  , 


f  =  -  M  log  ^^   +  f*  ,    P  =  -  Ma  log  —^   -  M  log 


r-1 


+  F* 


We  substitute  the  above  quantities  in  (2-3)  and  assume  that 

all  the  terms  involving  the  second  order  products  of  the  starred 

quantities  can  be  neglected.   The  linearized  equations  are: 

for   Q     <    X\     <    1      ,  -co<4<a)^ 

"I  =  -  Pf  OT)  '-    ^fPOTl   ^ 
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p*    -.  HAo*    : 
for     l<T|<r,  -oo<|<u-. 


(2.-0  U|   =   -    P|F^     :-F|P^    , 


pfe  =  -  ^^^v;;  +  ^*%)  -  p;; 


aU*   -I-  A*F^    +  AqP-  =   0    . 


P*  =  MAA*  . 
together  with  the  boundary  conditions 

f*(|,0)  =  0  ,  P*(|,,r)  -  0  , 

f*{4,l)  -  F*(|,l)  ,       p*(?a)  =  p*(?a) 

From   (2.'0   we   find  that   for     0  <  t^    <   1 

f|66    =   -   M"^A(MA    -    l)'-^(r-Ti)2ft         1-  2M"^A(MA    -    l)"^(r-ri)f; 


P|  ^  [(MA  -  i)''anr^(Pof|^  +  f|Po^)   ; 


and  for     1    <  t|   <   r    , 
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^ii  =  -  Ma(aMA  -  i)  (^"^^^^Inn  ■'"  m(J^  -  1)  (r-^)P|.n  > 

sub  ect   to 

f|(^,0)    =  0    ,  A|(^,r)    =  0    , 

f|(|,l)    -  P|(e;l)    >  p,f(l,.l)    -  ^|(?,1)    . 

The   solutions   of    (2.5)   which  satisfy  the   boundarj/-   conditions 

at     T)    =  0     and     r)   =  r     are 

s  s 

f|  =  a^   cos(v|    -;-  b)[(r-T])    "^    -   r2"^(r-ri)    2]    , 

p|  =:  a^    cos(v(i    :-  b)   ^(^^"]   -^ )    [-(M-^s^    :-  A)(r-Ti)    ^   •:- 

r2'''(M"^S2   •!■  A)(r-Ti)^l]    , 

^1 

Ff  =  a-   cos(v'4    H-  b)(r-ri) 

M~^S-^    :-  aA  S-|_  ;^1 


A|  =  -   a^   cos(v^    +  b)    rfa(^f.,^   .   1)    (r-n) 


,                               In.    ^(M>^    -    DM   ,.2.1/2 
where  (^  =  t^   [1    ' ^ ^^ —  v    J 

d  A 


s^    -  -  i    :-  CO   ,  Sg  =  -  I  -  CO   , 
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S^  =  I  [-  1  -:-  (1  +  ^("^^  -  ^^^^°  v2)l''2j  , 

and  a^  ,  a^  i  and  b  are  arbitrary  constants.   Application 
of  the  boundary  conditions  at  ^  =  1  yields 

-  (M~^s^  -:-  A)(r-1)  ^  +  r^'^ilf^s^  +   A)(r-1)  ^ 
..  [(r-l)'l  -  r2-(r-l)'2](M-ls,  -.  a.)   ^^^  "  ^^ 


1   ^"'  a(aM7v  -  1)  ' 

Now  v;e  define  the  critical  speed  £     as  the  limltinE  value 
of  A  when  v  — *-  0  .   Me  obtain  from  the  above  equation 

(1  -  M^)  =  {alU   -  i)(^  -  1)  , 


and 


0   =  ctr  ±  [(ar)^  -  !rr(a-l)  j-*-'^ 
2(a-l)M 


V/e  see  that  for  given  values  of  a  and  M  we  always  have 
tv70  critical  speeds.   It  will  be  shown  in  Section  5  that 
0<^M<1,    l<-iJ,M<co  if  l<a<««,   l<r<oo, 
where   i_   corresponds  to  the   "-"  sign  and  &    :-  .  the   "+" 
sign  in  the  expression  for  ^  ,   Futhermore .  if  we  rewrite 
the  quadratic  equation  for  Z     as 


am  .  ("-^^^^^^^  , 


it  is  seen  that  aM  -   1      is  alv;ays  positive  for  1  <  a  < 
and  I   <  r  <   03  .      Therefore   for  small  values  of   v  . 
AM  -  1  ^  0  and  aAM  -  1  >  0  . 
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Let  us  now  consider  the  case  a=l,   l<r<<».   First 
suppose  that  MA  ^^  1  -   We  obtain 

f|  =  P|  =  a^  cos(k^  -;-  b)[(r-ii)  "^  -  r^^(  r-n  )''^  ] 

which  is  unbounded  at  t)  =  r  as  seen  from  the  expressions 
for  S-,   and  Sp  .   Therefore;  the  linear  theory  falls,  other- 
v/ise  we  must  set  a-,  =  0  .   If  MX  ^  1  .  from  (2.5),  vie   obtain 

(r-Ti)f|^  -  f|  =  0  ; 

since  f|(|.0)  =  0  ,   f|  must  be  identically  equal  to  zero. 
The  above  results  confirm  what  we  already  established  in  Part 
I.   For  an  infinite  isothermal  layer;  either  the  linear  theory 
fails  or  the  solution  is  a  trivial  one. 

Finally  vie   shall  justify  the  consistency  of  our  lineariz- 
ing procedure.   Since   f?  ,  of  are  bounded  for  0  <  t)  <  1  ^ 
it  will  suffice  to  examine  whether  F|  .  A|  are  bounded  as 
T^  — ^  r~  .   The  terms  involv;ing  the  second  order  products  of  the 
starred  quantities,  vjhich  v;e  have  deleted  from  the  equations  for 
the  upper  layer  are : 

%   r\    ■      T]    %    '  5^     %    g  T]  Or) 


A^U*F*  .  A*F-  ,  A*U*F* 

0  T]    ^  T) 


We  see  from  the  solutions  for  F^  .  A^-  that  those  terms  are 
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of  the  order  (r-r))     or  (r-T])     and  tend  to  zero  as 

11  — ^^  r~     if  all'K   -  1  >  0  -  However,  for  small  values  of  v  , 

aMA  -  1  >  0  for  l<a<oo.   i<r<w..  and  the  linearizing 

procedure  is  then  justified. 


'i.  Nonlinear  Theorv.   Solitary  Wave  Solution 


We  assume  that  a  solitary  wave  moves  v;ith  a  speed  such 

that  X  =  -^  is  near  some  positive  value  &     which  is  to  be 
c 

determined  later.   Ue  write  (2.3)  in  the  form, 
for  0<'n<l_.   -oo<^<oo, 


(il.l) 


u.  =  -  P.f^  +  f^p^  , 

uffig  +  u.f.  -  (^-A)pf„  -  iof„  -  p„ 


puf„  =  1 


p   =   -   M(.g-A)p  -:-  lUp    . 

Similar  equations  v/ill  hold  for  the  upper  layer.   Now  let 

e  =  i  -  A  and  introduce  a  new  independent  variable  a  =  ^  >^ 

Then  ('41)  becomes, 

for     0<ri<l,      -co<CT<oo, 


u     =  -  p  f     -;-  f  p 


mo'    '".>-■   ni. 


:.:'--nT    .-,M       .-f-:--r;r    t~- 


.^^ 


f:% 


hi 


e(uf        ;■  u  f    )    =  epf     -   iof     -  p 


puf^^   -  1    , 


p   -   -   Mep    +  M^p    ; 
for     l<T)<r,      -co<0<co 


U     --p:^      +PP 
a  a  T]         a  Tj 


{k.2)  e(UP^cr   "^^  Va^    "  ^^^^t^   "  ^"^^^ti   "   '^^ 


AUP     =  a 


P   -  -   MeA   -;-  MM    , 

with  the  boundary  conditions, 

fCa.O)  =  0  ;  P(a.r)  =  0  , 

f(a,l)  -  P(a,l)  ,  p(a.l)  =  P(cr.l)  . 

It  is  assumed  that  all  the  quantities  in  the  new  independ- 
ent variables  can  be  expanded  in  integral  powers  of  e  ;  i.e. 


{H.3)  <J>(a.ri-e)  -  YZ     ^^'^J(^>'^)    > 

k  =  0 


vihere  <&  stands  for  p  ,  o  ,  u  ,  f  ,  P  .  A  ,  U  ,  and  F  . 
Substitution  of  ('4.3)  in  {'-i.2)   will  give  a  sequence  of  equations 


0     rr       I 


{■ 


-br 


•    c4M 


(P 


^  ari 
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and  boundary  conditions  which  these  equations  must  satisfy 
by  equating  the  coefficients  of  equal  poivers  of   e  .   The 
values  of  S,     and  also  the  solitary  wave  solution  will  be 
determined  by  solving  the  equations  for  the  successive  ap- 
proximations . 

The  equations  for  the  zero-th  order  approximation  are, 
for  0<T^<1,  -oo<a<c>o. 


0   -   -    ^PO^^OT)    -   POTI 


Po"o--oti  =  1  > 

('44)    for     l<ri<rj      -«)<a<co 


U        -_pp        +PP 


0   =   -   P^VOTI    -    % 


^O^O^On   =  °^   •' 


with  the  boundary  conditions 


"^.tr 


-c'l-: 


9rL1 


CO    s-    t)    >    oc 


.     O    - 


^3 


£^{0,0)    =  0  ,  ?q{o.t)    -  0  , 

fQ(a,l)  -  FQio.l)    .         Po{a,l)  :=  F^ia  ,1)    . 

We  may  assume  that  u^  -  U„  ;^  1  vjhlch  expresses  a  parallel 
flow  In  the  equilibrium  state.  The  solution  for  the  zero-th 
approximation  is 

Pq  =  i(r-Ti)  ,  Pq  =  ^(r-T))  , 

(-4.5)        .  -. 

Pq  -  (M)~^(r-r))  ,       Aq  =  (M)"-^(r-n)  , 

fQ  =   -   M  log  ^    ,  Pq  =   -   Ma   log  ^  -   M  log  ^   . 

The  equations   for  the  first  order  approximation  are 
for     0<'n<l,      -co<a<ooj 


^la   =  -   ^la'Ox}  '-  ^laPon   ' 


°  =  Po%   -   ^^Po^lTi  ■'   Pl%^    -  PlTi  ' 


"l   =  -   Po^lTi   -   Pl%   ' 


( 'I   6 )    for     l<T]<r,,      -oo<a<~. 


Uia  =  -   ^la^On  ^  ^la% 


0  =  ^VoTi  -  P^^Vm  ■^'^i^^  -  %  ' 


•.J    6   t; 


*  frJ  •        ^  pi-,    i'^^        p- 


,C'5>^)>00-.  , 


(TO     ..  ,. 


rrj;^  -  ^t^£'-- 
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«^1  -   -    ^O^lri  -  ^%  ' 


sub  ect  to  the  boundary  conditions 

f^(a,0)  -  0  ,  P^(a,r)  -  0  , 

r^{o,l)    =  F^(a,l)  ,     p(a,l)  -  P(a,l) 
From  (4.^)  we  have  by  elimination  of  u,  ,  p,  .  and  f,  , 

Piann  ^  °  ' 

and  then  f,   is  obtained  by 

The  solutions  for  o,_    and  i\ „     which  satisfy  the  condition 

'la         la  "^ 

at  T^  =^  0  are 


Pla  =   ^i^^^^M  -  ^^  ' 


r^^-.  riMai(a)^ii^5Jlil 


> 


where  a-!  (a)   is  an  arbitrary  function  of  a  .   Since  it  is 
assumed  that  tlie  flov;  reaches  the  state  of  equilibrium  at 
o   -  -  00         i.e.   a^(a)— 5-0  as  a  —^  -   oo   ^     we  have 


Pi  "  ^l^'^^^MI  ■  ^^ 


u  a 


0  ^    ( '•^  r- 


..-.;•/ 


D.t 


(r  -  ^)( 


(Mi    -    /)'r    ,_^  ._„    I- 


rr-1 


(o)].Bfl   •"l    - 


T)I 


n£    ji   9:>rit?      .    n     'io  rto 


5i       {  C  ) 


do 


((^ 


•1  - 
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f,  ^  -  r^  Ma,  (a)  ^^^  -  ^'^^)  , 
1  1       r-Tj    ' 

v;here  we  assvime  that  a,  (a)   is  not  identically  equal  to  zero. 

Similarly  the  solutions  for  A, _  and  F, _  which  satisfy 

la       la  "^ 

the  condition  at  t)  =  r  are 


A 


^^   =  Aj^(a){r-n) 


F^^  =  -  &  ^  m:^[(3)(-   1  -;-  ma)    , 


where  A' (a)   is  an  arbitrary  function  of  a  .   By  the  equi- 
librium condition  at  x  =  -  <»  ,  we  have 

A^  =  A^(a)(r-ri)  , 

F^  =  -  i"-^r-lAj_(a)(-  1  +  M.ga)  , 

where  we  assume  that  A.,  (a)  ^  0  .   Making  use  of  the  conditions 
at  the  interface  x\   =  \    ,     vje  obtain 

a^(a)(^  -  ^^  "  A^(a)(r-1)  , 
V\   7) 

r^   r4a^(a)  ~:^   =  "~'^   MA^(a)(aM  -  1)  . 

Since  A, (a)  ,   a,(a)  J  0  ,   it  follows  that 

(1  -  M^)  -  (^  -  l)(a^M  -  1)  , 

and  

n.  ON     s   ar  ±  i/(ar)   -  '!-r(a-l) 
^^  ^^     ^^  2(a-l)M • 


I  ■•'■-'^■ 


'Vi 


nJiBcrc! 


n     S' 


(  0 )  ^  A  -    (  X 


t    - 


X-'i 


.    (i 


m 


(C- 
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Thus,  we  have  two  values  of  i  which  confirm  the  two  critical 

speeds  obtained  from  the  linear  theory.  The  results  for  the 

first-order  approximation  are  summarized  as  below: 

(4.9) 

p^  =  a-^ia){^  -  r])    ,  Aj_  =  A^(a)(r-n)  , 

P3_  =  ma^{o){^   -  Ti)-(r-ri)  ,  P^  =  (MM^(a)  -  l)(r-T])  , 
f^  =  -  Z-^Ma^ia)    ^^   ;_^^^^  ,  P^  =  -rh]A^{o){am  -   1)  , 


u^^  =  -  Ma^(a)  ,  U^  =  -  MA^(a)  . 


In  order  to  determine  a-,  (a)  and  A  (a)  ,  we  must  proceed 
to  the  equations  for  the  second  order  approximation,  which  are 
found  as  follows: 
for  0<ri<l,   -oo<a<oo, 

"2a  =  -  (P2a^0Ti  +  Pla^ Iri  ^  "^  ^^2o^0r^   "^  ^IqPIti  ^  ' 


^laa   =    (Po^lTi    ■'  Pl^aq)    -    ^(P0^2T1   +  Pl^lri    "^  P2%  ^    "   P211    ' 


Po"2-OTi  +  Po"o^2-a  ■'-  p2"2^on  ■"  Po^i^iTi  +  Pi^o^in  +  Pi^i^on  ^  °  ' 


P2  =  ^^2  -  ^Pl   > 

(4.10) 

for     l<n<r,      -oo<a<w. 


U        =-PP        -:-PP        +PP 
"2a  ^2a^0Ti    '    ^2a^0Ti   ^  -^la^ln    ' 


Ls>i::iti.o  oWw    orii  r 


•  v;  \. 


t     (I    -   ?•■' 


V,  -■i)~(r 


f'.  -• 


-;     ,    ir>)^A     -■'      '■  ' 


:a1   ^' 


on    >     'D    >    00     - 


<  frx^o[^  ^ 


'  ps^i  -  (pQ-^s^  +  prx^  -^  fr<ro^^^  -  ^ 


'    '^   =    fro'^r'l^  "   nX-O^X^  -'   frx'^X^O^  "^   frO 


i    ^oV\  -  ^c' 


«  > 


,     >    X 


'    ,U't)X 


^^r"?      +  -^ 


OC 


U 
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together  v^ith  the  boundary  conditions 

f^{c,0)   =  0  ,  P2(CT,r)  =  0  , 

f^ia,l)    =  F2(a,l)  ,  PgCcT,!)  =  PgCa,!)  . 

For  the  lower  layer,  we  find  that 

(4.11) 

^2a  =  ^P0n^"^(^^P2a7i  -»■  ^'^S2a(^^^^  ^  ^^P20^On  "  Sl^^'^^^ 
where 


X  e3(j(^>^)  +  6iT^(c^.Tl)  > 


(4.12) 


g3(a,Ti)  =  -  p^f^^  +  uf 


\    'vi'\'^)£S      - 


(..' 


r '    .'i  > 


'i''. 


(r  '^l^ss- 


ri%i^  ■*• 


«  r^v>P  •"  r.c'*^^ 


•  x"-*-  n  I 
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Let  p  r\ 


Sll*io,r])    =  J      I  gi^(a,ri')dr!'d'n"     . 

V/e  have  from   (4.11) 

p2(j  =  M"-'-g|*(a,ri)   -  a^(cT)ri  +  b^(a)    , 

-  MifQ^a^(a)Ti  +  MKQ^b^(a)   +  r^g2^(a,Ti)   -  g^(a,Ti)]   , 

X'/here     aMcr)     and     b'(cr)     are  two  arbitrary  functions  of     a   . 
Since     f^   (a,0)    =  0   ,     by   (4.5)   it  is   obtained  that 

b^(a)  =  ^  a^ia)  -  ^  gf  (a,0)   -  ^  g|*(a,0) 


Hence, 


P2a  "=   ^m  '  ^^^2^^^   "^  M-\|*(a/o)   -  ^  gf(a,0)    , 

i'i      A/ 


(4.13) 


-  M-^gF(a,0) ^  gp„(a,0)   +  ^  g-.  (a,0)    , 


fg^  =  -   r^[-  rfe^(a)   -  M/.fQ^Tia'(cr)   +  rfQ^a^(a)   -i-  g|(c7,ri) 

+   ^fQ.,^g]l*ia,r])   +   rlg2^(a,r))   -  g^(a,T))   -  |  ^onS^CcJ^O) 

-  J^r^^s-f{a,o)  -  ^  fQ^g2^(cT,o)  +  r.  fo^g^(a,0)]  . 


.  f^•..-^ 


zo  aao: 


T  »■! 


',  T)J,\£. 


)     jTiO" 


.      ('V, 


V  ^A,  .I-/ 


(o.T7)|5,.ji  Jl  -  (rr,tj)j.vt  -  (p\t))^o:t    5>  +  (r 


.  Mo.o' 


(0. 


^9 


For  the  upper  layer  v/e  proceed  in  the  same  v/ay,  and  obtain 
that 

(4.1^) 

where 

(4.15) 

G,(a,n)  -^^1^%  -^Fla^m  ' 


G2(a,n)  =  U3_  -K  ml  -  m^^  +  f^^^  , 


G^(a,ri)  =  -  ^i^i„  +  ctU^  • 
Let 

G|f(a,Ti)  =  /   G^(a,Ti')dTi' 


p  n  r  n  " 
G|^(a,Ti)  =  J      I         gi^{a,Ti')d-a'dTT'  . 

It  is  obtained  from  (4.14)  that 

^2a   "  I^^G|*(a,ri)  _  A^(a)ri  +  B^(a)  . 

Since  A2^(a,r)  =  0  ,   it  follows  that 


-^Qrii 


I 


;^£ri;J   (4*.:. 

.    (c;  ^  -  ( 
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Ao„  =  A'{a)(r-ri)    +  M"-'-G**{a,Ti )    -   M"-^G**(a,r)    , 


'2a  2 

(^.16) 


Pga  =  Po^'-^[G*(a,T))   -  M^(a)   -i-  MiP^^A' (a)  (r-ii ) 


+  -^^^r^^'"^^   -   -^%^r(^'^^    '-  -^'^Gsa^^'"^^   -  ^1^^^^^^ 


The  matching  conditions  at   the   interface     ti   =  1 
then  give   tv/o   relations   betv/een     A' (a)      and     ai(a),    i.e, 

(4.17) 

M(^=-p5Y^)a^(0)    -  M(-   1   -:-  M^a)A^(a)    =C  J^o)    , 
where 

Ci(a)    =  M'^G|J*(a,l)    -   M"^G|*(a,r)    -  M~^g|*(a,l) 


+  M'^gj*(a,0)    +  -^  g/f(a,0)    +  — ^  gpa^^'O^ 


-|-    gn  (C7,0) 


Mi 


C  2^^^    =  -   S/t^^^l)    +  ^  g|(cJ,0)    -  ii£y  g|-*(a,l) 


+  frr  sr  ^^>o^  -  ^'Sa^^'^^  -^  TfTUT  e2a(^'0) 


-1-  g3_(a,l)    +  ^  gi(a,0)    I-  G|(a,l)   +  f§  G|[Ma,l) 


^  G^(a,r)    +  r^G2^{a,l)    -   Q^^o ,\)    . 


(I 


{0,D'. 


iL,v)^- 


.  .••<■' 


^^■^•Qs^  rc^^  ^^-^-t.;- ■"■ 


'i'- V 


.-•'  'p 


{s:,x}]'^p  ^      (i,d)^o 
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From  the  results  of  (4.7)  ,    the  following  relation 


h  -   ^  r-l       ^_1^ 

.„!   -  m  -   M(aM  -  1)  qAo) 
^     r-l  "2 


must  hold  In  order  to  ensure  a  consistent  system  of  equations 
for  aA(a)  and  A'(0)  .   Let  S,     assume  the  values  given  by 
(4.8)  .   Then  we  have 

(4.18)    M(a^M  -  l)6\(a)  =  (r-l)C2(a)  . 

By  some   complicated  but   straightf orv/ard  computations  and  re- 
arrangements  of   the   terms,   we   finally  obtain  the  equation 


where 


mQaj^''(a)   +  m^a-[(a)a-^(a)   -^  m^a.}^{a)   =0   , 


.  ^^'^^   -   f )    [2(l.r)r(l   -  aM)los  ^ 

0  (my  ^'^ 


1   +  r  •;-  m{l  -  2ar)]    , 


(4.19) 
m 


T    =  -^^  \  4r(r-l)log  -^   [(a-l)(M)^   -  a(M)^  +   (^M)^] 

^      {my  *^  ^'^ 

-   2r(a-l)(M)^   +   (-   2r^a  +  6ra  -h  a  -   r  •:-  1)(M)^ 


+  (ar^-   5ar  -i-  4r^-  lOr  +  4)(M)^+   (4ar^-   3ar  -   3r^ 


+  12r)M  -    (6r^-    3r)j    , 


M^rr2 
m^  = 


2  ■  IW^IM  "    (2-Ki)] 


•70  I 


T0~'>[ 


J.-  :■- 


X   •:•   1   - 


f  v^  r  .  ii  ) 


•i<L 


-  tl»5;    -"'rai' 


u^..    /        - 


.  I 
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Let  us  assimie  that  m^  ,  m.,  ,   and  nip  are  not  equal  to  zero 
for  given  values  of  r  ,  a  ,  and  since  we  are  only  concerned 
with  the  solitary  wave  solution,  the  following  conditions 

a^(-  co)  =  aj^'(-  ex,)  =  0  ,    a^(0)  -  0 

are  imposed.  The  solution  for  a, (a)  subject  to  the  above 
conditions  is 

(4.20)    a^(a)  =  -  -^   sech"^  | 

1        \ 

Then  by    (4.7)   v/e  have 


nig 


r 


(,.,a,    .^,„,  ..«_^^3ec.^|L^. 


Nov;  suppose  that  the  successive  approximations  so  far 

obtained  do  furnish  a  sufficiently  accurate  representation 

of  a  solitary  wave  in  a  compressible  medium  of  two  isothermal 

layers,  we  obtain,  in  terms  of  the  independent  variables 

X  and  7  , 

for     0<r|<l,      -oo<fi<oo, 

o                  r,                3mp            p     -,      fm   {-K-&) 
p  ^  ^^gH^(r-Ti)-  ^^C2(i-A)M^(^  -   r,)  -^  sech^  '^ J'\ ' 


p  ^  M-^   Pi(r--q)-   Pi(i-A)(fe  -  r^ )  -^  sech^  -^/-^ 


0 


u  ~  C  -1-  MC(^-A)  — ^  sech^  ^  Z-^- 

m^  2H     /      mQ 


7  ^  .  ►m^iog  x^  .  r^MH,  nO_z_M)  (,.,,2^  ^eoh^  ^  '"^^^ 


2; 

mQ        ' 


^^ 


'    0",-.C' 


>     V  t-f: 


c::»j.<J  !>^' 


V  ^'' 


o«^ 


( '^. 


IT- 


^  IT  '  rio  r 


[A-t) 


jr-'i 
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(1   -  M^)3m^     ^  /m^(A-^)  ^     „        /ni^(A-i;) 


V  ^  Ci'^M ^  -JI-(i-A)/_2 gg^j^2  _x_     y_2_ ^ 


m,  r-n  /      m^  2H,     /      m. 


for     l^T^isr,      -co<a<oo, 
(^.22) 


'-^2,,.^,  ^,,,  II   -    ^    5^2  .2      X       I^^P^^ 


—  1  -^  -   1   3^0  o      ,      /rn-(  A-^) 

A  ^  M-%(r-,)-   ^^(...)(r-n)  ^%n--^  -^^'  #^/-  ' 


mo 


^  -   1    3mo  o     -      /m^TXrn 


U-C-HMC(£.A)M__Z^3ech2^j! 


niQ 


F  ^  -   MH^a  log  -^  -   Ml-I^log  ^  +  r^m^    ^\v.l)^''    ^^"^^   ^ 


jiriglA-i!) 
X       ~ 


niQ 


r-l  m^  niQ  2H^    /      m^ 


r/ 


3^       im^tA-i!) 
X  tanh 


2H,    J      niQ 


Where  ^ 

p,  ,  gH-, 

Ti   =  -2—   ,      R     ^  -^  [exp(-Sl_)  -   1]    ,      A  =  -^   , 
PjCH^                        gPl                Pl/Pl  ^ 


r  =  1   +   [exp    (— S^)    -   1]"-^    ,      M  =  r  -   1    . 

Pl/Pl 


7  can  also  be  expressed  In  terms  of  the  vertical  distance 
^  at  X  =  -  00  .   For  0  <  r]   <  1    , 


r  -  "• 


s.::-;i::f5i.&    Isolc^isv    srlJ 


<    X  2: 
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0  Pi        Pi 


and  for  1  5  fl  <  r 

^.  c      gp^  h         p  C 
-  [exp(-   " 

^^1  h 


7  =  -^  [exp(-^^— )  -  1]  -1-  /   p^u^dy 
^        Pi  ' 


PnC       gp,h  gp 

-4-[exp(— ^)  -  exp(-  — ^  (C-h))]  . 


Oy 


Pi  P] 


5.  Properties  of  the  Solitary  Waves 
Near  the  Tx-io   Critical  Speeds 

As  mentioned  in  the  Introduction ^  it  can  be  seen  from  the 
complicated  coefficients  of  the  solitary  wave  equation  that 
a  complete  analytical  approach  to  the  study  of  the  solution 
for  the  solitary  waves  is  almost  impossible.   However,  v;e 
shall  try  to  explore  analytically  the  properties  of  the  solu- 
tion as  much  as  we  can  with  the  help  of  numerical  calculations, 
To  begin  v/ith  v;e  first  indicate  some  interesting  properties 
of  the  critical  speed  Z    .      Prom  (4.8)  we  let 

p   M  -  cxr  -:-  [(ccr)^  -  4r(a-l)]^/'^ 

(5.1) 

0  ir^  ^   ar  -  [(ar)^  -  ^r(a-l)]^^^ 
2(a-l) 


-'d.. 


asJ. 
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vjhere 


T 
2 
00  >  a  =  TfT-  >  1 

^1 


00 


>  M  =  r  -  1  =  [exp(-^^^ — )  -  1]"-^  >  0 

Pq/Po 


vie  claim  that  oo  >  ^j/i  >  i  and  1  >  i_M  >  0  .   This  can  be 
shown  as  follows:   rewrite  the  expression  (ar)  -  4r(a-l) 
in  the  square  root  of  (5.1)  as   (ar  -  2)  +  4(r-l)  ;  hence, 

(ar)^  >   (ar)^  -  4r(a-l)  >  0  . 


This  shows  that  ^_M  >  0  .  Now  let 

z  =  \  -  m  . 

The  quadratic  equation  for  I 


becomes 


(a-l)(M)^  -  ar(M)  -i-  r  -  0 


(a-l)(l-2)^  -  ard-z)  +  r  =  0  , 


i.e.      {a-\)z^   +  z[-  2(a-l)  -i-  ar]  +  {a-l)(l-r)  =  0  . 

However,   a  >  1  ,   r  >  1  ,   the  two  roots  of  z  must  be  of 
opposite  signs,  i.e. 

z   =1-^,M<0,   z   =1-^M>0, 

since   iM,  >  iM_  .   This  completes  the  proof. 

Next  we  shall  investigate  certain  limiting  values  of  ^ 


-:-   I  o>  - 


.   0  =   (i-I){.C-x>)   +   ['IX) 


lo   9c'  C'Eurn     ■? 


,   C  <  M^i,    -   I   =   _i:      , 


■.'J 


■-1  J. '-■■.". 


lOO'X'i   ocll 


J      lo  ast.'iflp,' 
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Assume  that  M  is  positive  finite.   Then, 
as  a  — >'  1  '  , 


^  M  _>  ^-  CO  ,  £   n  -^  1^    ; 


as  a  — >-  00  , 


i_i_M  -t-  r"^  ,      i_M  -^  O"*"  . 


Nov;  assume  that  co  >  a  >  1  .   We  have, 

a. 

as  r  -^  1 '  , 

£  M  -^  _i_'  ,    if  a  <  2  , 
+     a-1 


1^  ,   if  a  >  2  ; 


^  M  ->  1"  ,    if  a  <  2  , 


1 


a-1 
as  r  — >*  00  , 


,   if  a  >  2  ; 


With  the  knowledge  of  the  properties  of  the  critical 
speeds,  we  proceed  to  investigate  the  coefficients  of  the 
solitary  v/ave  equation  given  in  (4.19).   It  is  seen  that  the 
solution  of  equation  (4.19)  with  any  one  of  m^  ,  m,   and 
mp  vanishing  v;ill  be  identically  equal  to  zero  if  the  con- 
ditions a,  (-  <»)  =  a-[(-  co)  =  a-I'(-  «)  are  imposed.  There- 
fore, vje  must  determine  along  which  curves  in  the  domain 


3iA?     ''O     ......  .        - 


►T'.ii       r'"    <    .jf"' 
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a  >  1  ,   r  >  1   those  coefficients  will  vanish.   V/e  take  up 
nip  first  which  is  given  by 

2 

Let  mp  =  0     and  assume   that     l    <  r  <  co     and     1   <  a  <  oo    . 
Then  vje  have 


a+2  • 
Substitution  of  the  above  equation  for  iM  in 


(5.2)     {a-l){m)^  -  ar{m)    +  r  =  0 


y^^l^^        Ma-1) 


r  = 


.2   ,, 


a^  -  4 

The  portion  of  the  locus  of  the  above  equation  in  the  domain 
a  >  1  ,   r  >  1   is  shovm  in  Fig.  5;  and  &  =   &_     along  the  curve, 
below  v/hich  mp  <  0  ,  and  above  which  mp  >  0  .   For  the  values 
of  ^ ,  ,  mp  is  alv/ays  negative. 

The  coefficient  of  aj^''(a)   given  by  (4.19)  is 


mQ  = 


^'^   ^^   ~J^'^^    [2(l-r)r(l  -  aM)log  -^  +  ^   + 

+  ^M(l  -  2ar) ]  . 

11(1    -    2M) 
For     oo>r>l,      03>a>lj     — ^^ — /  0     and 


(M)^ 


m^^  =   0     implies 


2(l-r)r(l   -  a.m)log  ^  +  1   -h  r  +  M(l   -   2ar)   =   0    . 


-jcU  r  i^vCs 


s        * 


q      .   0  <  ^in 


2i   (01,4*)    -'Id  ao.  Z-)      (c)  "  j-S 


1    >-    I  -f 


bac     0 


»v  r 


im) 


^    X   <   jD  •<    ••» 


(ir>2   -    i:::1  -;■  ■": 


58 
In  combination  viith    (5.2)   v;e   find  that 

[2r(r-l)log  -IL.  _   1]   ±   i  [_   2r(r-l)log  -^  -)-  if  -  4r(r-l)[ 

M  =  ^^ '■ ^^ 

2[2r(r-l)los  -^  -   2r  -;-  1  ] 


r-1 


2r(r-l)log  ^  -   2r  -f  1]/ 


If    .1.1 


{m)^  -  r           2r(l-r)log  ^  -;-  1   +  r  +  M 
a  = = —  =  

(M)^  -   r(M)  2Mr[(l-r)log  ^  +  1    ] 


Of  course,  vje  may  eliminate  ^M  and  obtain  an  expression 

for  a  in  terms  of  r  ,      but  it  is  much  more  involved  than 

the  relation  betvjeen  Jil'l     and  r  .   There  are  two  curves  found 

in  the  domain  r  >  1  ,  a  >  1  along  v/hich  mQ  =  0  (Fig.  5)« 

One  corresponds  to  £^    ,   and  the  other,  S,      .      Let  us  denote 

the  former  by  £_-cur\'e     and  the  latter  by  <g, -curve  .   Then 

it  is  shovm  that  m^  >  0  belov;  the  £   -curve  and  to  the  left 

of  the   >g, -curve  ,      and  m^  <  0  above  the  2   -curve  and  to 
+  u  - 

the  right  of  the  -2  -curve. 

Finally'',  v;e  are  confronting  the  coefficient  of  a.(cr)a-,(a)  , 
m^  ,   which  can  be  rev;ritten  in  the  form, 
for  co>r>l,  co>a>l, 

m^  =   (r-l)[4r(r-l)log  ^  [a((M)2  -  (M)  +  1  -  {M)^] 

+  r^a[-  2(M)  +  1   -v   4(M)"^]-:-  ar[-  2{m)^  -v   6(M)  -  5  -  3(^M)"^] 

+  r^[4  -  3(M)"-^  -  6(M)"^]-i-  r[2(iJM)^  -  (M)  -  10  +  12(M)"-^ 


,-2 


+  5(M)"^  ]  +  am   +  (M  -:-  h)f 


7 


•IJ.ro- 


.     (D).c(tj) 


t   X  <  o 


-  I  +  im ) 


(^"(it:^)c  -  p  -  m.)<^  •; 


.[-. 


((•r  )G.'i  +  or  *  (MS.)  - 
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where  J  by  (5.2) 

(M)2  -  r 


a 


{UDim  -  2)  • 


We  find  by  numerical  calculations  that  m,  >  0  for  Jl  =  & , 
and  there  is  a  i_-curve  in  a  >  1  ,  r  >  1  below  which 
m,  >  0  and  above  vjhich  m,  <  0  . 

Let  us  denote  the  six  subdomalns  in  y  >  1    ,     a  >  1    , 
divided  by  the  three   i_-curves  by  D_-I  ,   D  -II  ,'" ,   D  -VI 
and  the  tv;o  subdomalns  by  the  only  i  -curve  by  D,-I  ,  D,-II 

"r  "r  T 

(Fig.  5).   We  collect  our  results  as  below: 


m-[_        nig 


Domain 

mQ 

D_-I 

+ 

D_-II 

^- 

D_-III 

- 

D__-IV 

- 

D  -V 

- 

D_-IV 

- 

D  -I 

+ 

+ 

D  -II 

+ 

+ 

+        + 
+ 

+ 
+ 


In  order  to  determine  v/hether  the  solitary  vjave  is  of  elevation 
or  depression  type,  \ie   must  go  back  to  (4.22).   It  is  seen  from 
the  expressions  for  f  and  F  that  for  given  M  and  r  the 
wave  amplitude  increases  as  t)  increases  and  reaches  a  maximum 


(1) 

Wa' 

ve  Type 

+ 

D 

- 

D 

+ 

D 

- 

D 

- 

D 

+ 

D 

+ 

E 
E 

(1)   The  sign  of  H-'h     must  be  such  that  the  expression 
mp 
-  {H-'K)   — ^  is  alv;ays  positive. 


c 


a 


VT        r 


TT 


ri'^  '  -•    'io    of    9'  .  ■.    ofl;^ 
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at  T^  =  1  ,   then  stays  at  constant  value  for  1  <  "H  <  r  . 
The  wave  type  is  determined  by  the  sign  of  the  expression 
(1  -  m){£--K)m^ 


m^ 


V/e  denote  by  E  the  wave  of  elevation 


and  by  D  the  I'jave  of  depression,  and  the  last  column  in  the 
above  table  indicates  the  wave  type  in  each  subdomain. 

Novi  let  us  consider  the  limiting  case  that  M  is  posi- 
tive finite  but  a  -:-  1"*"  .   Then 

(i^M)  -^  CO  ,     (^_M)  ->  1"  . 

The  latter  is  the  only  case  we  need  to  consider. 
As   (i_M)  ->■   1"  , 

mo->  0  , 

m3_  ->  (r-l)[-  2r^  +  2r  -h  6] 

m^  -^  '   r(r-l)   , 

hence  there  exists  no  solitary  v/ave  solution.   This  limiting 
case,  in  fact,  shows  that  as  the  temperature  of  the  upper 
layer  and  lov/er  layer  tend  to  the  same  value  the  solitary  wave 
will  disappear  as  vje  have  observed  before. 


ym:'''j  .'>■ 


.0      lO     9V 


.' ods 


..     -.-W    pp" 


•<—    V  1"»    * 
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Part  III.  Compressible  Media  of  Infinite  Depth 
with  Non-Uniform  Velocity  Distribution  at  Equilibrliun 

1.  Introduction 

In  Part  I  and  Part  II  \-ie   are  concerned  with  the  situation 
that  a  solitary  v/ave  Is  generated  by  some  disturbance  in  a 
compressible  medium  initially  at  rest.   In  reality  this  is 
usually  not  the  case.   In  the  equilibrium  state  the  atmosphere 
may  move  with  non-uniform  velocity  distribution.   Tlie  vjork 
done  here  is  to  investigate  the  solitary  waves  in  a  compressible 
medium  initially  v;lth  arbitrary  velocity  profile.   For  sim- 
plicity, we  only  consider  a  medium  of  Infinite  depth  at  con- 
stant temperature.   Nevertheless,  the  method  employed  here 
can  be  extended  to  cases  of  multi-layer  polytroplc  or  Isothermal 
compressible  media  of  finite  or  infinite  depth  without  much 
difficulty. 

The  formulation  of  the  problem  is  presented  in  Section  2. 
In  Section  3  the  linear  theory  predicts  a  critical  speed  for 
an  arbitrary  equilibrium  velocity  distribution.   We  treat  the 
problem  by  means  of  the  nonlinear  theory  in  Section  4  as  vje 
did  in  Part  I  and  Part  II.   The  solitary  vjave  solution  is 
computed  based  on  a  velocity  profile  of  exponential  growth. 
It  Is  confirmed  that  there  exists  no  solitary  solution  as  the 
velocity  profile  tends  to  a  uniform  one. 


r 
J. 


'.VJ 


vf  Ibl'c'. 


Icit   ;«■;.' i: bam 


-  lu 
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2.  Formulation  of  the  Problem 

We  consider  a  compressible  medium  at  constant  temperature, 
which  fills  up  the  whole  upper  half  space.   The  medium  is  sup- 
ported bjr  a  rigid  plane  bottom  and  the  pressure  p  is  asstmied 
to  be  zero  at  infinity.   There  are  no  geometric  constraints. 
In  the  equilibrium  state  the  velocity  profile  as  a  function 
of  the  vertical  distance  only  is  assumed  to  be  given.  A  cross 
section  of  the  medimi  is  the  upper  half  plane  (Pig.  ^).   It 
is  assumed  that  a  v;ave  of  permanent  type  moving  to  the  left 
with  constant  velocity  c  has  been  generated  in  the  medium 
due  to  some  disturbance.   V/e  choose  a  coordinate  system  m.oving 
with  the  wave  such  that  the  x-axis  coincides  with  the  bottom 
and  the  y-axis  passes  through  the  crest  or  the  trough  of  the 
vjave  and  is  positive  upward  (Fig.  k) .     V/ith  respect  to  the 
coordinate  system  the  v:ave  is  stationary  and  the  velocity 
of  the  medium  at  Infinity  in  the  equilibrium  state  is  denoted 

by  u^(y)  • 

The  governing  equations  are 

a(p^)  .,  $i(pv)  ^  0 

d::   '    dy      ' 

~  Su  ,  ~  c)u     1  ^p 
dx     oy     ~  ox 


'V  oV     ~  OV 

u  -^r-  ■:-  V 


1  Si 


^     '^   ~  oy  ' 


"Sx  "■'  '  "Sy   ~     "^  °y 


•^^^'^     '■>  I" 


9'J  • 
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=  p 

Po 

v;here  p(x,y)   is  the  density,  u(x,y)  ,      v(x,y)  are  the 
horizontal  and  vertical  velocity  components,   p(x,y)   is  the 
pressure,  g  is  the  gravitational  constant,  and  p^  ,  Pq 
are  respectively  the  reference  pressure  and  density  ivhich  may 
be  taken  as  the  ground  level  values  in  the  equilibrium  state. 
Denote  by  ^(x,y)   the  stream  function  such  that 

fy  ^   pu  ,    fy^-   pu  . 

The  mass  ilvzz   across  any  vertical  plane  from  y  =  0  to  y  =  «> 
per  unit  breadth  is  given  by 

p   00 

^(x,co)  -  ^(x,0)  =  /    P^u^cly  =  7^ 

^  0 

where  we  may  set  ^(x,0)  =  0  . 

Let 

^{x,y)  -  7 

where  0<y<co,   0<7<7   ,  and  assuine  that  for  any 
value  of  7  in  0  <  7  <  7   there  exists  a  unique  solution 
of  the  above  equation  for  y  such  that 

y  =  f (x,7)  . 

Hereafter  the  bar  notation  is  used  to  denote  a  quantity  as  a 
function  of  x  and  7  .   If  we  choose  f  ,  p  ,  o  ,   and  u 
as  dependent  variables  then  v;e  have,  as  obtained  in  Part  I, 
for  0  <  7  <  7   , 


\^ri»  'loT   ^iQf'.i    e;:j;t:;".s  i  ri:.;      ,      r  >   ^  >   0 


i.-'/pci;;  B  V 


Oo 


i>si{.7   rloaa     1^     r' 


»    ( 


,       '/»      >        '      >     nil     — 
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"x  -  -   ^Px  ■'-  ^xPy  ' 


(2.1) 


P  u  f^  -  1  , 


Po   Po 


subject  to  the  boundary  conditions 

f(x,0)  =  0  ,    p(x,~)  -  0 


In  order  to  make  v.       and  p   appear  explicitly  in  the  above 
equations  it  v;ill  be  more  convenient  to  use  the  vertical  dis- 
tance C  of  the  stream  lines  at  x  =  -  oo  as  an  independent 
variable.   The  relation  between  7  and  C  Is  given  by 

(2.2)     7   =J        Pco^cc^y  > 
0 


and  it  follows  that 
(23) 


A 

dct)    _   1      ^cp 
G„ 


vjhere 


^  -v      /  «  \~ 


G=PjOuja    , 


CO       ^  CO 


and  the  hat  notation  denotes  a  function  of  x,  C  •   However, 
if  the  transfomation  defined  by  (2.2),  (2.3)  Is  always  pos- 

A  A 

sible,  we  must  impose  a  condition  on  G  ,  i.e.   G  is  non- 


zero. Since  0  (C)   is  alv/ays  positive  for  0  <  ^  <  »  , 

'  CO   ^ 


,  c 


.'i3vrO'.7oH      .    ■;.    ,x     'io 


«   >    ^   >   0      •ro'V    -v-JLv^i.Scq    os^svrffi   ai      ( 
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we  assiome  that 


for     0  <  ^   <  CO   , 

B^   using  jc,  ^  as  independent  variables  (2.1)  becomes, 
for  0<^<oo,   -oo<x<<», 

A  A        A  A      A  A 

G  u ,  ==  -  f  „p   +  f  p„  , 

AAA       A  A  A  A 


^co(^  ^\x  -^  Vx^  =  -  PSf^  -  P^  , 


(2.4) 


AAA 

A 

p    U    f^    = 

G 

A              A 

£_  _  P 

<V                   "V 

^ 

Po        Pn 

subject  to  the  boimdary  conditions 


f(x,0)  =  0  ,     p(x,oo)  =  0 


In  order  to  non~dimensionali2e  the  above  equations  vje  introduce 
the  following  dimensionless  variables 


A 

A 

1  = 

X 

h   ' 

V 

^ 

V 
c    * 

A 

u   = 

U 
C      ' 

A 

Tl     = 

^ 

P 

= 

PqC 

p  = 

f    = 

A 

f 

h    ' 

X 

= 

2    ■' 
c 

Goo 

A 

Geo 

PqC 

Poo 

0 

•    00 

""  -..     •» 

Po 

"co 

^0, 

~    c       ' 

t     h'-- 


.     0     =     (c--. 


.   0 


.    '^--^  u 


^  71 


-  I 


=^     0 


r,  • 


=-  A 


C>1 


u 


O 


ee 


0 
where     h  =  -^^—  .      Then  v/e  have,   from   (2.^), 


0 

for     0<ri<oo,        -co<a<co, 


00  I;  -q^^  '^^ri 


(2.5) 


G^(uf^«    -I-  u,f,)    =   -   pAf^   -   p^    , 


p   u  f„    =  G 

'  T)  CO 


P    =    Ap     , 

with  the  boundary  conditions 

f(^,0)  =  0  ,    p(^,co)  =  0  . 

These  equations  look  quite  similar  to  those  in  Part  I  and 
Part  II.  We  expect  that  the  same  procedure  given  before  may 
as  well  be  applied  to  the  present  problem  vjithout  introducing 
much  difficulty.  Hov/ever,  since  u  (ri)  is  arbitrary,  it  is 
anticipated  that  the  coefficients  of  the  solitary  wave  solu- 
tion will  be  much  more  Involved  than  before. 


3.  Linear  Tlieory.   Critical  Speed 

VJe  first  recall  that  in  the  equilibrium  state  the  density 
distribution  for  an  isothermal  medium  of  infinite  depth  is 
of  exponential  decay,  i.e.. 


(  n )    SJ      r. 


-0 


6? 


In  terms  of  p   and  rj  ,  we  have 
(3.1)     Pco(^)  =  exp(-  Ti)  . 

The  steady  state  solution  for  a  parallel  flow  with  velocity- 
profile  Uq  =  ^^('"i)   is  found  from  (2.5)  and  (5.1)  as  follows: 

Pq  =  ^Pco(n)  ' 

(3.2) 

fo  =  -n  . 

Suppose  that  a  small  disturbance  is  superposed  on  the  parallel 
flow  and  vje  may  write 

u  =  u^  +  u*  ,      p  =  Xp^  +  p^'^ 

p  =  p^  -I-  p*  ,      f  =  n  +  f*  . 

If  we  substitute  the  above  quantities  in  (2.5)  and  assume 
that  any  term  vjhich  contains  second  order  products  of  the 
starred  variables  can  be  neglected,  v;e  have  the  following 
linearized  equations : 

G  u  f,1?,.  =  -  X(p  f*  +  p*)  -  p*  , 


^r 


»1  -}-    r:  =■-   1 


5Vf     , 


<     f  oo^'^' 


.    fq    -    (♦q  + 
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(3.3)     P*u  -1-  0  u*  -:-  G  f*  =  0  , 

subject  to  the  boundary  conditions 

f*(^,0)  -  0  ,    p*(^,~)  =  0  . 

It  is  obtained  from  (3«3)  that 


G  u  Tf-.-.    =  A(u^  -   A)'-^?  u^f,1:^„ 

00    CO    l;L;t^  co  i^oo    co    (;T|T1 


2 
'  u     -   A  ' 

CO 


and 


o*  =   (u^  -   A)"-^[Ap   £f  -    0  u^f;?-    ]    . 

r  ^     CO  ■■     "^co    t  •  00    CO    t  n 


lie  first  try  to  find  a  solution  for  ff  .   Let 

s 

f|  =  H(UP(t))  , 

s 

and  from  (3.4)  v/e  obtain 


is 


Tne  solution  for  PI  is 

H  =  A  cos(v'-  -V   B) 


(.  i . 


A    -       iJ 


.    0    - 


f.i      .   f'l     -ro"^. 


.    0 


-^ae 


.■^t.ik^)'"'^,^   ^  p(S^')'-U  -^..)  -1  .  ^p? 


■■(:  -^<) '-,-.%  - .  i„(f,.-.o^  ^- 
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where  a  and  b  are  tvjo  arbitrary  constants.   Since  the 
critical  speed  i     Is  defined  as 

(1) 

i  =  lim  A(v)  , 
V  ->  0 

the  follovjlng  asymptotic  method  is  suggested  to  find  the  value 

of  Z     vjhile  a  general  discussion  of  the  solution  of  P  vjill 

be  given  in  the  Appendix.   Let  us  assume  that,  for  small  values 

of  V  , 

p     h 
A  =  i  -I-  V  A^  +  V  Ag  +  •  •  •  i 

The  equation  for  P^(r) )  is  found  as 


2    ,^-l,.,2^   ,  /„  .,2^-l/^  ,,2 

(3.4) 


%^  +  [-  (<  -  ^)'-'(<)-n  -'-  (Peo<)    (Pc»<)J% 


subject  to 

F(0)  -  0  , 

lim  (u2  -  ^r^ihj^   -  9 Aw   =  °  • 

r\   — *-  CO 

The  general  solution  for  Pq  is  obtained  as  follows: 
Fq  -  C[.  1  -;-  ILe"-    j        e'"^ 'u^^dV  )dri '  ]  +  De'l  . 

n 


(1)   For  definition  of  the  critical  speed,  c.f.  [1]. 


fiil-J  tzi'.LZ      .^jii-^.jaoo  -iii^'icJicliE 
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The  condition  at  n  =  ^  is  no  more  than  a  boundedness  condition 
for  P  and  requires  D  =  0  .   Since  ^q(^)   =  0  ^  ^^   we  assume 
that  the  motion  is  other  than  a  parallel  flow,  i.e.   C  ^'  0  ; 
\ie   obtain  the  critical  speed 

£  =   [T^   e"''"''u~2{-T  )diV  ]'^  . 
0 

It  is  interesting  to  consider  the  case  u^  =  const  . 
(1)   u^^  7^  A  . 

Prom  ( 5 • 5 )  we  ha  ve 

p    -  p   =  -  v^-K'-^iu^    -    A)P  . 
r)T)    'n  '  c:.    ' 

The   solution  for     P     vjhich  satisfies      P(0)    =0     is 

m,  r,  mpT] 

F  =  c(e   ^      -   e   "^    ) 

v;here      c      is  an  arbitrary  constant  and 


m,    =  I    [1    -  Jl    -   ^v2(u2   -   A)X-^    ]    , 


mp  =  ^   [1   -:-   Jl   -   ^^v2(u2  -  ajA"^    ]    . 


VJe  may  use  the  condition  p?(ii,")  =  0  to  determine  A  and 

2  2 

consider  the  cases  u  >  A  and  u  <  A  .   Hoivever  P  is 

CO  :j 

ali^ays  unbounded  as  t\  —>   co  except  that   c  =  0  .   Hence  for 
u  7^  A  linear  theory  fails  or  the  solution  for  P  is  a 
trivial  one . 
(2)   u^,  =  A  . 

From  the  eouation  for  op     in  (3-5)  we  have 


ijyxjx.. . 


■J         .0.1         YOJ 


.  [ 


:o'^   sjfj^H      .0-0     :ja:l:i  3. 


rr 
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P„  -  P  =  0  . 

It  is  seen  that  F  =  0  if  P(0)  =  0  .   The  above  results, 
in  fact,  are  exactly  v;hat  we  obtained  from  the  linearized 
equations  for  the  case  ti  =  1  in  Part  I,  except  that  a  differ- 
ent independent  variable  is  used  in  each  case. 


4.  Nonlinear  Theory.   Solitary  Wave  Solution 

It  is  assumed  that  a  solitary  wave  moves  with  a  speed 

such  that  A  =  -^  is  near  some  positive  value  i   ,  which 
c 

is  to  be  determined  later.   The  equations  (2.5)  cai^  t>e  writ- 
ten as, 

for  0    <    T\     <    CO     J        -c»<^<co, 

G  U/.  =  -  f„P/.  +  fgp^  , 

G  (uf,.,.  -1-  u.f.)  =  (i-A)pf„  -  ^pf„  -  p„  , 
(4.1) 

p  =  -  (i-A)p  +  Up    , 

subject  to  the  boundary  conditions 

f(4,0)  =   0  ,    p(e,«')  =  0  . 

I7e  let 

e  =  £-\ 

and  Introduce  a  new  variable 


XT 


JO    9V. 


o    n^z\'   .;■'">■ 


(r.ii) 


t     l^ 
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Then  (4.1)  becomes, 

for  0<'n<oo,   -oo<a<oo, 


(4.2) 


eG  (uf  „  +  u^f  )  =  £pf  -  ^pf  -  p„  , 


puf  =  G  , 

"^   T)      CO  ' 


p  =  -  ep  -h  £p    , 

together  v;ith  the  boundary  conditions 

f(a,0)  =  0  ,    p(a,oo)  =:  0  . 

As  before,  suppose  that  all  the  dependent  variables  can  be 
expanded  in  integral  powers  of  e  ,  i.e. 


(4.3)     ^{o,r],e)   =  Z_  e^  fv(a,ri) 

k  =  0     ^^ 


where   <!>  stands  for  p  ,  p  ,  u  and  f  .   Substitution  of 
(4.5)  in  (4.2)  v/ill  yield  a  sequence  of  equations  and  bound- 
ary conditions  for  the  successive  approximations  by  letting 
the  coefficients  of  like  pov;ers  of  e  be  "equal. 

The  eqviations  for  the  zero-th  order  approximation  are, 
for  0  <  T\   <  CO   ,      -  00  <  a  <  CO   ^ 


0  =  -  ^Po%  -  Pon  ' 


(4.4) 


io    CIO:  ......      .3        .     1       .     — 

■bnuofi 
r  vd   ail'.'/  ■ 

,  ....    ,   .      30'       3       10     _. 


^    eo    >    X)    >    <»    "       t« 


'    [70*^*00"  '■    ;t:j'oO^ 
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PO  =  ^9o   > 

with 

fQ(a,0)  =  0  ,    p(a,~)  =  0  . 

Let  us  assume  that  u^-,  =  u  {r\)     v;hlch  represents  a  steady 
parallel  flovj.   The  solution  for  the  zero-th  order  approxi- 
mation is  found  as : 

Pq  -  hj-^) 

Po  =  Pcc(^) 

fo  =  ^  • 

The  equations  for  the  first  order  approximation  are, 
for     0<Ti<co,        -oo<a<t»j 


CO  la  Or)^la         la^Qn    ' 


(^.5) 


(4.6) 


0  =  Po%  -  ^(Po^m  '-  Pi%)  -  Pin  > 


PiVoTi  +  Po^i%  +  Po^o^in  =  °  ' 


Pt    =    -Pi    -    Po    > 


'1        ^0 
vjith  the  bovuidary  conditions 


f^(a,0)    =  0    ,  P^(cr,oD)    =  0 


It  is  obtained  from  (4,6)  that 


,    '-.' 


/  ^ 


n  =. 


(^.^) 


0     ii-ol 


-  .• ; 


.li  rictc 


.    0   -    («"tD)j.q  t    0 


7^ 

Integration  o±    the   above   equation  and  making  use   of     p-.  (cJ^")    =  0 
jj-lelds 


r  "  P,  (''V  ) 
P.  -  =  a  •  ( a )  /       -^ dTi  ■ 


where     a' (a)      is   an  arbitrary  function  of     a    .      Since   the 

flovj  is  assumed  to   reach   the   equilibrium  state  at     x  =  -   "   , 
i.e.      a(a)— s^O     as     o  —*-  -  co   ^     v/e   obtain 

pcopJV) 

p,    =  a  (a)/        -^ dri'    , 

where  we  assvme  that  a (a)   is  not  identically  equal  to  zero, 
It  follows  from  (4.6)  and  the  solution  for  p,   that 

=  -  r\-\^{a)[.   Z   ^^^  •;-  /   >-^— -  dV]  . 

By  the  equilibrium  condition  at  x  =  -  -^  ^   vje  have 


f-,    =  -   ^-^0  -^a(a)   -   &  ^p.^  •;-  /       -^ 

1  ■    '•-3  '    C3  /  d  I        ,    \ 


Since     f^(a,0)    =  0   ,     and     a(a)  £  0,      it  follows   that 
-   1   +  i  -^/       -5 dri>    -   0   , 

'o    ^^^('^'^ 


"■rlw 


f  ' 


El".' J'    "^woXXol    i^i       ,0    ._ 
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CO   p    (r) '  )  -1 


r,  CO     Oil-; 

i.e.  i  ^    [         -^ dTi'] 

^0      ^^(^'^ 


The   value  for     £     confirms   the   critical  speed  we  have  found 
by  the   linear  theory. 

In  siimmary,    the   solutions   for  the  first  order  approximation 
are 


(4.7) 


Pt    =  a(a)/       — 5 d'q  •    , 


1  ^     "•  U  00 

CO 


f,    =  a(a)[r2  -     1      /       -i^^ dn '  ]    . 

1  ^Pco^'  v2(^.>) 


To  determine  the  function  a  (a)  \ie   must  go  one  step 
further.   The  equations  for  the  second  order  approximation 
are,  for  o  <  'o  <  "  ,   -  oo  <  a  <  co  ^ 

Geo^^oo^^aa^  -    (Po^n  ■■'-  Pl-%^  -  ^(Po^Vo  ■'■-  Pl^lr,  "^  P2%^  "  "^2^    ' 

Po^2^0-o  +  P0^'0^^2-n  ■•■  !^2'-0^^0n  "^  Pq^I^p.  ■■  Pl^^O^^Ti.  -'-  Pl^^l^^OT.  =  0  ' 
(4.8) 

Pg  =  ■^?2  "  Pi  ' 


c   if  j. .Co ire D   9ficr 


{■r.5. 


nb  — 


/r    jl 


T.ii) 


[  'rrl) 


]iD) 


c?9j»;  ero 


.T9:J;>b  cT 


Oft 


nv 


--I) 


0 


<     U9^    -    ^f.o'^^.'-   '     frl^L^   -'• 


0    - 


rro^i"!^  ■'■  m^o'^r;  ^  ux'r^o 


•  rUr^Q     -J- 


t ;« .  i- 1 
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subject  to  the  bo\.indary  conditions 


It  is  obtained  from  (^,^.8)  that 
2    ^    ?-  2 


^^-P2-naN  -  ^  "ooPoaTi   =  Si,(a,n) 


^20   =  Pop"^    f '"^"^P2Tia   -'■  P2a--0vi    '  ^'^'^2a  "   ^i  ]    ^ 


v;here 

-:-  K%a'^    '■  P0r,TiP0Ti~^fS2a  "    ^^i  1   -  -^Sp^    , 
gl(c7/.l)   =  -  p^^f^^   -;-  f^^p^^    , 


cop      (i-j'  ) 


Now  from  the  condition  f,^(a,0)  =  0  and  -2  =  [/  -^5 d'o '  ] 

we  obtain  the  follov/ing  condition 

(4.10)    u2(0)p2,^^^(a,0)  V  P2^{a,0)[r'   g     diV  J""^  -:-  Sga^^'O) 
-  gT(a,0)[/   -5 dTi>]--'  =  0. 


The   equation  for     Pp        can  also  be   rev;ritten  as 


(4.11)     (^^)   '^e^(a.n) 


CO  /  'CO 

11 


'  I'i 


^"^'s- 


iv; 


'['[it 


\  . 


ra^xc/xfjf^- 


O    3W 


(0.1') 


.    0    ^ 


0 


-    ad   o?.  1 


Dii^^ 


.   (fft^) 
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Multiplying  both  sides  of  the  above  equation  by   /    —^ 6.r\ ' 

and  then  integrating  v.'ith  respect  to  'q   from  0  to  oo  ^   v;e 
have  ,  .  < 

0      CO   • 

0      T)     CO   ' 


It  follov;s   from   (^.10)    that 


/         (/        -| d-nM-i 


0       ^      CO 
r     ?co(n') 

=   /       -p clTTg      (a,0)   -  g.  (a,0)    . 

CO        ' 

By  some  lengthy  but  straightforward  calculations,  from  (^.9) 
and  (4.11)  vie   finally  reach  the  equation 


mQa'''(a)  -;-  m^a '  (a)a(CT)  -i-m^a'Ca)  =0 


where 


(1)   In  deriving  this  equation ;  v;e  must  assume  the  condition 

lim   ( /   —^ drj '  )  — >■   =  0  . 

n  — !*  CO  ■'    u  (n '  )        ^co 

T]  CO   ' 

p 

However,  from  the  equation  (4.11)  if  u^  y^  0  as  r\   — *>  co  , 

we  find  that  gu(a,T])  =  0(6"')   and  P2ar\   "^  ^  ^^  ^1  — »'  0  , 
and  the  above  condition  follows  immediately. 


no.. 


-  \ 


0  - 


:io. 


■L'^ieflw 


/T     as 
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"^0  =  -/''  P>(^)f^"^cj^^f  -  ^"MF(n)]^dTi 


^  0 


-  F(0][r^G{0)u^{0)    -    ^"^u^(0)F(0)]  , 


0  ^00 

-  (2  -;-  p"^)F^(r|)  +  4p"V^F(-n)  -  o"^  (p"^f2(ti  ) 

-  2u;2F(Ti))^]dTi  -  r^F^co)  -  i-V2(o)  , 


mg  =/   P(ti)[-  4u;2  .  P(ri)(p^  -;-  p'J-)]&^-]   -  r^F(O)  , 


F(ri)  =  /   -g dT]'  . 


Suppose  that  none  of  itIq  ,  m,   and  nu  are  equal  to  zero, 
and  v/e  also  Impose  the  conditions 

a'  (-  od)  ^  a'  '{-  oo)  r.  0  ,    a'  (0)  =  0  , 

then  the  solution  for     a (a)      is 


a(a)   =  -  -^  sech^  |     -  -^ 
m^  2  J     itIq 


Ass\ime  that  the  successive  approximations  up  to  the  first 
order  give  a  sufficiently  accurate  representation  of  a 


nb(i   - 


,.T 


f ;  r 


OH 


=    ,nf 


c    r„      r. 


.  (o)-^-^-'i  .  (o)H'-"^  -  .^fc[„((f-)'^ 


{0)'d~*S.    -    f't.r  (^"q  -;-      v,;)(p-).^ 


f;v    -  i  t 


0 


-  ;  JJ  i'l  •  pill 


(.a    -)'  -.3    =     (co    -)  'S 


a  J:      'r.)G     'lo': 


•riJ  ndrl;^ 


0  S 


0'-'     i-  - 


o.) 


1 ' 


Sail' 


£  lo  noJtii 


■  a-itoLl'^jje. 
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solitary  vjave.  1/e  have;,  in  terms  of  the  Independent  variables 
X  and  ^  ,      Tor  0<^<c3,   _oo<x-oo. 


;?m. 


1  ^0 


m. 


1  -i      0 


^  ^  r  ^"S  _  n-l^-l. 


f  =  hTi  -  h(i-7v)[^"''  -  r-^p;^F('n)]  -^  sech^  f  j-  ^  (i-A) 


0  -^  [r2  .  rlp:>(,)]  ^^^4^  (.  ^  (i.A))V2  3eeh2  X 


"7 


h 


m^ 


v/here  ti  =  -^  ,  0 

'   h    ■  c 


Po 


exp(-  r\)    ,      h  =  — ^^  ,   A 


c 


5.  An  Example 

In  the  follov;ing  v;e  shall  give  a  concrete  example  based 
upon  a  special  velocity  profile 

u  =  1  -  k  exp  (-  T])  ,     0  <  k  <  1  . 


Tiie  tv/o  extreme   cases  ^     k  —>■  0'     and     k  — **  l"    ,   v;ill  be 


ysesi. 


.  I- 


A   -   L 
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discussed  at  the  end  of  the  section.   By  repeated  integrations, 
we  find  that 
m. 


v/here 


^  -  ^Ic^  -   5k  -   5 
6(l-k)^ 


m     =  1 ^   [1^5  .;.  5j^^  _   27k5  +  32k^  -   25k  -^  2]    , 

^       6k(l-k)^ 


2k  -   1 
2        2(l-k)^ 

i  =    (1-k)    , 


hIq  <  0    ,  for     0  <  k  <  1    ; 


m,    >  0    ,  for     0  <  k  <  kQ    , 

=  0   ,  for     k  -  kQ  =   O.OS    , 

<  0    ,  for     ICq  <  k  <  1    ; 


rog  >   0   ,  for     1/2  <  k  <  1    , 

=  0    ,  for     k  =  1/2    , 

<  0    ,  for     0  <  k  <  1/2    . 

The   expression  for  f     is   found  as : 

^^  ^  C   -  h(^-A)[l   -  (l-k)(l   -  ke-^)-l]   ^^^^f^^Y  ^^   sech2  f 


J3(2k  -  l)(4k2  -h  5k  +  5)('2-A) 
where 


N(k)    =  k^  -;-  5k^   -   27k^  -i-  52k^  -   25k  -:-  2   . 


.    [S 


::&c  -:    -:ITS  - 


'-A-D. 


X    - 


.    0    < 


5\I  >  >1  >  0     •10'^ 


ai     1     'lo 


"  ''"'■  — nw    '    ^   ^^' 
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For  a  given  k  ,   the  v/ave  amplitude  Increases  as  -q   increases 
and  finally  reaches  a  finite  value  as  ri  — »-  «  .   The  v/ave 

nip 

type  depends  upon  the  sign  of  -  (i3-A)  — ^  . 

Let  us  divide  the  open  interval  0  <  k  <  1  into  three 
open  intervals,   I^  =   (0,kQ)  ,   I^  =  (kQ,l/2)  ,  and 
I-^  ==  (1/2,1)  .   IJe  list  our  results  as  follov;s: 


I  , 

niQ   , 

m-j^    , 

nig   , 

^1 

- 

+ 

+ 

^2 

- 

- 

-1- 

I, 

^^ 

_ 

i-X  ,  Wave  type 

+         E 
+         D 


V/hen     k  =  k^   ,      k  =  1/2   ,      there  exists  no  solitary  v;ave   solu- 

A 

tion.   It  is  also  seen  from  the  expression  for  f  that  as 
k  -^  0  ,  f -j^-s-  0  ,  and  as  Ic  -^  1  , 

f-^-:-  +  (^-A)  |j  sec^  g  ^42(;2-A)  ,   for  t]  >  0  , 

A 

f^-5^  0  ,   for  r|  =  0  . 

The  former  confirms  that  the  solitary  wave  solution  disappears 
for  an  isothermal  layer  of  infinite  depth  as  the  velocity 
distribution  tends  to  a  uniform  one.   Tlie  latter,  however, 

A 

shov/s  a  very  queer  situation.   The  limiting  solution  of  f  as 
k  — >-  1  possesses  a  discontinuity  at  r\   ^  0    .      This  is  phys- 
ically Impossible  and  the  theory  fails.   In  fact,  in  this 
case  m^  ,  m,   and  m^     tend  to  +  co  or  -  oo  as  k  — >-  1~  . 


-.3.1  i         .     ■•    -^ —    i' 


ij.'ia 


A  .. 


^     0    <    (v       to'l        t    \a'-''  '' 


a'i 


3  ay  J 


0    v:-- 


a.!:f'..J    ai    ^: 
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Appendix  I.  Solution  of  the  Linearized  Equations 

in  Part  I. 

V/e  shall  investigate  the  solution  of  the  equation 
(I.l)        ^  1 

_  1 
n 
subject  to  the  boundary  conditions 
P(0)  =   0  , 

1-i  i 

lim   (l-ri)[(l-Ti)   "  -  ^r^[(l-^)\  -  AF]  =  0  , 
Ti  ->  1 

In  order  to  facilitate  our  discussion  we  let 

nil 

2  =  (l-n)  "   , 

and  (I.l)  becomes 
(1.2) 


.n-lx2, 


n  n 


The  two  independent  solutions  of  (1.2)  for  v  =  0  are  found  as 


(1)   Hereafter  we  use  the  same  capital  letter  F  to  denote 
a  function  of  z  . 


rtoi::^£5;jp''/   ■"- '. 


.C 


-.    i         p." 


r-  i 


.  Q  -  o .. .  .  r 


.?.fo/.-t:.  :nMJj&d  eii;/  oJ   ■Joei.dua 


.•f.L 


(i 


(   r     -\      ;-,rtf5 


f :";.!) 


8!?  bnuo'i   S-TJE     0  =--  V     tcl   (S.I)   lo  BfioiSL-i..^  -  ... 
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P^  =  1  -  Az  , 


Now  we  construct  two  Integral  equations 


F^  =  Fj_  -  v^j        K(a,z)F^(a)da 


(1.3) 


where 


Fg  =  Fg  -  v^r   K(0,z)F2(a)dcT  , 


K,         ^        F.^(z)F3(a)  -  Fg(z)F-^(a)      ^        ,  . 
K(a,z)  = a      (^^)  (a  -  ^) 

P^(a)F2(a)  -  F^(a)F^(a) 

1 

.:  .  ^^_^^^   [(1  .  ^z)   -  (1  -  Aa)(z-^a)"-^]  . 


By  the  theory  of  Volterra  Integral  equations  it  is  easily 

shown  that  since  F^  ,  F^  ,  K(a,z)   are  bounded  for 

0<'n-^<Ti<l,  '<]   <  o  <  1    ,      P,  (z)   and  Pp(z)   are  uniquely 

determined  by  (II. :5)  and  also  satisfy  (1.2).   From  (1.3)  we 

see  that   z  =  — r-  is  an  apparent  singularity  of  (1.2)  for 
nA 

0  <  ■=r-  <  1    .      It  is  recalled  that  the  equation  for  p;^  is 
nA  ~  '^ 

given  by  -,  -j^ 

p|  =  G{a(i-Ti)[(i-n)   "  -  HA^"^f^^"^^\  -  ^^^ 

n-1 
=  G((-)z  "  (z  -  ^)"^[-  ^  Fz  -  AF]  . 


(1.4) 


i  -■ ; 


('0 


.n 


)      D 


{0).^{0)~'^    •.    (O)^-?! 


,)?!    , 


I  .1 


Id   (S.I) 


[■T/^  - 


r^ 


.    ['^A  -   s^ 


X-r. 


f  t' .  r 
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for  V  =  0  then  it  must  be  so  for  v  ^  0   .      First  we  note  that 


VJe  claim  that  if  z  =  — r-  is  an  apparent  singularity  of  p^ 

n  A  c, 


(1.5) 


Let 


n 
n-1  „     ,„      ^~   n-1,     1  % 


F  =  c^F^  -I-  c^Fg  , 

and  substitute  F  in  (1.4).   It  is  seen  that  z  =  — r-  is  an 

nA 

apparent  singularity  of  pf  as   v  =  0  for  0  <  — r-  <  1  . 

For  V  ^  0  ,  we  rewrite  K(a,ri)  as 

1_ 

K(a/a)  =  [F-j_(z)F2(a)  -  Fg(z  )F.j_(a )  ]  (^_^)^  a     ^'^  , 

and  note  that 

K(z,z)  -  0  . 

Therefore,  from  (1.3)  v/e  have 

\,   -   F,^  -  v2/  [F,^(z)F2(a)  -  F2,(z)F,(a)]  j^  a^'^ 

z 

F^(a)da  ,  1=1,2. 


Let 


F  =  c^P^  -:-  CgFg  , 


where  c,(l-A)  +   Cp  =  0  since  F(l)  =0  . 
Then  from  (1.4)  it  is  obtained  that 


•J  ^ 


(-^ 


off 


.    0  ^    (i, 
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n-l 
p?  =  G(|)z  ^    (z   -  ir-)-^   X 

S  HA 


n-l    Ti      /_  \        -v  -n    /_  \  \     .    _    /      n-l 


-   v^/^f-  ^K^  Pi,(-)    -   XF^(z))P2(a)    -    (-  ii^  Pg.Cz) 
z 

]■ 


-   AF2(zl)P3^(a)]   ^  ~x  cj""^    (c^F^(a)    +  c^F^{a))da 


It   follov;s   from    (1.5)    that      z   =  -^     for     0  <  -^  <  1  is  an 

•^  nA  HA 

apparent  singularity  of  pfr     for  v  ^^  0  . 

s 

In  the  neighborhood  of  z  =  0  ,  v/e  apply  the  method  of 
Probenius .   Put 


m  =  0 


Substitution  of  F  in  (1.2)  yields 


aQ[(^)c(c-l)  +  c]  ^  0  . 


This  gives 

1 


c  =  0  and  c  =  - 


n-l 


and  v;e  ali^ays  obtain  two  solutions ,  one  of  vjhich  is  bounded 
and  the  other  of  vjhich  may  contain  a  logrithmic  terra  and  is 
alvjays  iinbounded.   It  is  easily  shovm  that  only  the  bounded 
solution  satisfies  the  boundary  condition  as  r|  — ?-  1~  given 
in  (I.l)  .   The  recurrence  formula  for  the  series  expansion 


<.>i 


L'. 


^uh- 


-    0 


a-:  : 
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of  this  bounded  solution  is  given  l)y 

2  2 

n       n        n  A 


J. 


or  m  s  0  ,  and  a_„  =  a_,  =  0  ,  and 


m  =  0 


(1.6)      =  aQ[l  -  Az  +  v^P(z,k^,A)]  . 


It  is  also  easily  shown  that  the  series  expansion  converges 
for   |z|  <^   . 

Finally,  vje  shall  verify  v/hether  the  linearizing  procedure 
is  consistent.   The  terms  v;e  have  neglected  from  the  full 
equations  (2.17)  are 

f*P|  ,  f|P*  ,  u-q,   ,   u|f|  ,  p-f-  ,  u-p-f^^  ,   u*Pof*  , 


p*f*  ,  p*u*f*  . 


Since  in  the  neighborhood  of  ri  =  1  , 


T)    z  d  0     \   I ' 


P,^  =  O(l-Tl)   , 

s 


■f  » 


,  r 


i. 

fl 
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u*  -  0(1)  , 

f?  -  0(1)  , 

It  Is  easily  verified  that  all  the  terms  listed  above  are 
bovmded  at  ■'"1  =  1  and  our  linearizing  procedure  is  then 
justified. 
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Appendix  II.  Soliitlon  of  the  Linearized  Equations 

In  Part  III. 

In  the  Append!::  vie   shall  consider  the  equation 


F_   -r   i-    iv.^  -   A)"^(u2)„   -:-   (p  u2)-l(p  u^)^]F 
■q  rj  CO  '    00  r|  '  co  cd         *  >  co  co  t)      t] 


(II. 1) 


■'■'  (^l  -  ^)"'(^v^)"'  ^(^^-f)-/  -  -  ■'^^^l '  ^^)^"'  p 


subject  to  the  boundary  conditions 
P(0)  -  0  , 

lim   (u^  -  A)"-^[Ao  F  -  0  u^P„  ]  =  0 

CO  -co      '  CO  CO  n 

T)  — >-  00  ' 

as  given  in   (3.5)  vjhere     u     ^  constant    .      The  tv/o  independent 
solutions  of    (1)   as      v  =  0     are  found  as 


-1       -ir~  Pco(^') 
^1  =  -  ^      -;-P^7     -277;^^'   ' 

■^  U     (T)  '  ) 

(II. 2)  ^ 

2  •   CO 


00 


vie  remark  that  the  equation  (1)  when  v  =  0  is  exactly  of 

the  same  form  as  the  one  for  f-,   v/e  might  have  derived  from 

the  set  of  equations  for  the  first-order  approximations  in 

the  non-linear  theory;  hov/ever,  these  tvjo  solutions  can  be 


TTSD-caiioo 


c?nsfc.-i-:  owd   ?^if]' 


(I. II) 


^^  ^-.    .  x- 


io  vj.cfo&xs  ci     0  -  V     aoriw   (I)  .••: 


•.ivi 
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easily  obtained  from  the  solution  of  p,   which  is  governed 
by  a  much  simpler  equation  than  f-j^  .   From  the  knoijledge  oJ 
the  reduced  equation  of  (1)  we  can  construct  two  integral 
equations 


or -a 


P3^(ri)  -  P^(ii)  -  v^/   K{a,r])F^{o)da    , 


(11.3) 


PgCn)  =  PgCn)  -  v^jT   K(a,ri)F2(a)dCT  , 

0 

X'jhe  re 

FAri)F^ia)    -   FAc)F   {a)        ,   p 
K(a,n)  =  -^ ^ ^ ^ A  ^(u2  -  A) 


(II. 4) 


¥^{a)F^[a)   -  F^(a)P3^(a) 


=  -  A-V(a)  4-  A-^Ja)p;l(Ti)u2(a)  h-  p;^(Ti)u2(a) 


T)     CO   ' 

For  any     r)      >   0     such  that      0<r;<TiQ<oo,      0<a<r], 

P,  (t))  ,  Fp(r))   and  K{a,Ti)   are  bounded,  by  the  theory  of 

Volterra  integral  equations,   F-|^(-q)   and  Fp(n)  are  uniquely 

determined  by  (II. 5)  and  satisfy  (II. 1).  As  seen  from  (II. 5) 

and  (II.4),  u^  =  A  is  an  apparent  singularity  of  (II.l) 

for  any  bounded  ri  .   We  also  recall  that  the  equation  for 

P*  is 

9t   =  H(e)(u2  -  -KV^i-K^J  -   9^^l\)    . 


■voji    ■■  •        •  ■  -yj 


Q       r.- 


03  CC    * 


)  ,'^. 


•■.II.' 


t    ;"   >    ^   i   •^      .-■''•    -y^   i    f   i   0      Jsrfcf   rl0JJ3      0   -' 
'to   'iioadd    9,i,t   Yo'    , 
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It  follows  from  (II. 3)  and  (II. 4)  that 


p|  =  H(^)  ';c^[y   -^ dri-  +  v'^pjn)j        K^(a/q)F^(a)da] 


L  -^  «.(n') 


.,.rv. ,...,. .,.,], 


-  C2[l  -  V   pj^)j        K^{a,ri)F^io)do] 

0 
Inhere 


-1   /_^_-l,.,^,.2,,^  ,   -1..,.  2.,,/"^  P-^^  ) 


K  (a,n)  -  ^■^oJa)p;^(1^)u^(a)  -;-  p;^(-q  )u^(a)  /   -^-—  dq  ■  , 

-^  ^     u  (n  '  ) 


and  C-,  ,  Cp  are  two  arbitrary  constants;,  and  by  P(0)  =  0  , 

T    p  "  p  (n' ) 
c,  (-  A"-^  +  /   -^ dri')  -I-  c„  =  0  . 


This  shows  that  u  =  A  is  also  an  apparent  singularity  of 
pf  for  any  bounded  rj  . 

At  infinity,  the  equation  (II.l)  may  not  possess  a  bounded 
solution  for  F  corresponding  to  a  given  bounded  u^(ii)  . 
In  order  to  facilitate  our  discussion  \ie   introduce  the  var- 
iables : 

G(t])  ^   e''1p(Ti)  , 


and  (II.l)  becomes 


nli'. 


•r.f)^-/^V 


I-. 
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V/e  assume  that 


and  let 


n  =^  0 


00 


G(z;    =     /^ a    2 

n   =  0     " 


Since  the  method  of  Frobenius  is  x-;ell-known,  in  vjhat  follovjs 
V7e  only  Indicate  the  general  results  v/ithout  going  into  details 
There  are  several  cases  xve  need  to  consider: 

(1)  bo  =  X  ,  b^  .  b^  =  •••  =  b^^_^  =  0  ,  b^  ^  0  . 

In  this  case  the  two  roots  of  the  indicial   equation 
are   c  =  0  and  c  =  m  -:-  1  ,   One  solution  of  P  is  bounded 
and  of  0(e~'  ')   and  the  other  contains  a  logrithmic  term 
and  is  of  0(e'')  .   It  is  shown  that  only  the  bounded  solu- 
tion satisfies  the  boundary  condition  at  infinity. 

(2)  bQ  ^  A  ,   bQ  =  -  b^_^  ==0,  \  ^  0    ,      m>0. 

(a)   bp^  ;t'  0  .   The  indicial  equation  is  found  as 
c2  .  c  -1-  v^{^   -  1)  =  0  , 

and  the  two  roots  are 

C  ±  =  I  [1  ±  03]  , 


^-_k.,.2\^ 


We  find  two  independent  solutions  for  F  free  of  logarithms 


-.1 
.  "I 


,     [iO    t 

A 


or:!o' ill ';y,o  I   'io  eoil     '^     'rol   a; 
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(|  -  f  )n 

if  o)  is  not  an  integer.   One  solution  is  of  0(e      )  , 

the  other,  of  0(e      )   .   If  b^  <  A  ,  then  to  >  1  ,   and 
only  one  of  the  solutions  is  bounded  and  satisfies  the  condition 
at  infinity.   If  b^  <  A  ,  co  <  1  and  both  solutions  are  un- 
bounded. 

For  m  >  0  ,   the  indicial  equation  becomes 


c^  4-  (m-l)c  -  v^  =  0 


and 


c  ±  =  ^  [-  (m-1)  ±  J'(m-1)2  ■;-  ^v^  ]  . 


For  small   v  ,  both  solutions  of  P  are  unbounded. 

Finally  v;e  note  that  all  the  terms  ive  have  neglected  in 
the  linearizing  procedure,  i.e. 


u*p^f*  ,   o*f "U   ,   p*u*f*  are  bomided  if  there  exists  a 

'^  0    T\  •  i^co'    ^        COT) 

bounded  solution  for  P  at  -q  =  co  .   This  justifies  our 
linearizing  procedure  given  in  §3,  Part  III. 


nl 
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Figure  f 
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Figure  2 
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D^-n 
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Figure  3 
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Figure  4 
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